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ABSTRACT 

The  Calculation  of  Mid-Latitude  Meridional 
Neutral  Winds  from  an  Improved  Servo 
Model  Approximation  Method 

by 

Theodore  R.  Ballard,  Master  of  Science 
Utah  State  University,  1991 

Major  Professor:  Dr.  Kent  L.  Miller 
Department:  Physics 

The  theoretical  servo  model  of  Rishbeth,  Ganguly  and 
Walker  provided  the  impetus  for  a  physically  more  practical 
version  of  a  servo  model,  called  the  USU  servo  model.  The 
USU  servo  model  was  developed  to  automatically  calculate, 
with  a  personal  computer,  meridional  neutral  winds  of  the 
upper  thermosphere  for  any  mid-latitude  location  from  F2- 
peak  layer  heights. 

F2-peak  layer  heights  were  input  into  the  USU  servo 
model  from  three  different  sources  for  calculation  of 
meridional  neutral  winds.  The  International  Reference 
Ionosphere  provided  mid-latitude,  geomagnet ically  quiet, 
monthly  averages  of  F2-layer  heights  for  different  solar 
cycle  conditions.  The  International  Reference  Ionosphere 
also  provided  a  means  of  comparison  with  ionosonde  F2 -peak 
layer  height  deduced  winds.  F2-layer  heights  from  differ¬ 
ent  mid-latitude  ionsondes  were  determined  by  using  scaled 
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ionosonde  critical  frequencies  from  an  empirical  equation 
developed  by  Bradley  and  Dudeney.  Millstone  Hill  incoher¬ 
ent  scatter  radar  provided  a  means  to  compare  the  USU  servo 
model  with  other  neutral  wind  calculation  techniques, 
including  radar-derived  winds. 

Results  of  this  investigation  indicate  that  the  USU 
servo  model  is  a  valuable  tool  in  the  calculation  of  mid¬ 
latitude  meridional  neutral  winds  and  compares  well  with 
radar-deduced  wind  data  during  quiet  to  moderate  geomagnet¬ 
ic  activity.  When  combined  with  F2 -layer  heights  produced 
by  the  International  Reference  Ionosphere,  the  USU  servo 
model  can  become  a  meridional  neutral  wind  prediction  tool, 
computing  solar  cycle  and  latitudinal  variations  in  neutral 
wind  circulation  patterns.  (117  pages) 


CHAPTER  I 


INTRODUCTION 


1.1  Overview 

This  study  involves  mid-latitude  meridional  neutral 
winds  derived  from  ionospheric  F2-layer  heights.  General¬ 
ly,  the  F2-layer  peak  (hmF2)  occurs  between  200  and  600  km 
[Hargreaves,  1979].  At  these  altitudes,  neutral  wind 
velocities  are  normally  constant  with  height  due  to  strong 
viscous  drag  effects  [Richmond,  1983].  Therefore,  winds 
derived  at  hmF2  are  regarded  as  constant  with  height  and 
approximate  the  neutral  wind  of  the  upper  thermosphere. 

Meridional  neutral  winds  are  calculated  from  an  im¬ 
proved  version  of  the  F-region  servo  model,  known  in  this 
study  as  the  USU  servo  model.  The  servo  model  has  been 
successfully  used  in  previous  studies  and  can  reproduce  the 
behavior  of  the  mid-latitude  F2-layer  quite  well  [Ganguly 
et  al.,  1980].  The  servo  method  provides  a  computational 
advantage  by  describing  the  F-region  through  physically 
simple  processes  at  the  F2-peak.  The  most  serious  assump¬ 
tion  is  that  the  model  deals  only  with  the  F2-peak,  so 
calculated  values  are  assumed  to  apply  throughout  the  upper 
thermosphere  [Rishbeth  et  al.,  1978].  Historically,  the 
servo  model  has  compared  well  with  other  methods  to  deter¬ 
mine  neutral  winds  [Buonsanto,  1986,  1990]  and  [Bjonsanto 
et  al.,  1989,  1990].  Buonsanto  et  al.  [1989,  1990]  found 
that  when  the  servo  model  is  "tuned"  from  deduced-radar 
winds,  servo  model  wind  speeds  are  closer  to  radar  wind 
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speeds  throughout  the  period  of  study.  Tuning  allows  for  an 
empirical  adjustment  of  a  balance  between  chemical  and 
diffusive  loss  of  plasma  at  the  F2-peak.  The  USU  servo 
model  is  divorced  from  empirical  tuning  and  is  physically 
more  realistic,  while  maintaining  the  servo  advantages  of 
speed  and  simplicity.  This  study  explains  improvements  in 
the  USU  servo  model  and  shows  how  the  model  may  be  applied 
to  calculate  the  meridional  neutral  wind  from  a  variety  of 
sources  explained  in  Section  1.3.  Calculated  neutral  wind 
patterns  are  also  qualitatively  related  to  momentum  sources 
during  geomagnetically  quiet  and  active  conditions. 

Global  studies  of  neutral  wind  patterns  that  respond 
to  changing  momentum  sources  can  improve  our  understanding 
of  thermospheric  dynamics.  However,  systematic  observa¬ 
tions  of  neutral  winds  are  limited  to  a  sparse  network  of 
Fabry-Perot  interferometers  (FPI)  [Hernandez  and  Roble, 
1976],  incoherent  scatter  radars  [Salah  and  Holt,  1974; 
Wickwar  et  al.,  1984],  and  satellite  measurements  of  the 
Doppler  shift  of  the  airglow  spectra  [Herraro  et  al., 

1988].  It  is  logical  that  we  can  improve  our  understanding 
of  thermospheric  neutral  wind  circulation  patterns  by  using 
a  method  which  can  expand  our  global  coverage. 

1.2  Problem  Statement 

Two  methods  of  deriving  the  meridional  component  of 
the  neutral  wind  from  hmF2  values  have  been  successfully 
used.  One  method  [Miller  et  al.,  1986,  1987,  1989]  uses 
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the  Field  Line  Interhemispheric-Plasma  (FLIP)  ionospheric 
model  developed  by  Richards  and  Torr  [1985].  The  other 
method  [Buonsanto,  1986,  1990]  and  [Buonsanto  et  al.,  1989, 
1990]  uses  the  servo  model  of  Rishbeth  et  al.  [1978].  The 
servo  model  has  a  computational  advantage,  while  the  FLIP 
model  is  much  more  comprehensive.  The  servo  model  is  more 
readily  adaptable  to  parametric  investigations,  that  is  the 
variation  of  neutral  winds  due  to  gas  composition  and 
temperature  changes . 

The  two  models  essentially  use  the  same  method  to 
derive  neutral  winds  and  therefore  should  produce  about  the 
same  results.  A  recent  paper  by  Buonsanto  et  al.  [1989] 
made  a  comparison  of  the  two  models  using  Millstone  Hill 
incoherent  scatter  (IS)  radar  data.  The  servo  model  pro¬ 
duced  winds  consistently  more  equatorward,  especially 
during  geomagnetically  active  periods,  than  the  Miller  et 
al.  [1986]  method.  During  this  comparison  the  servo  model 
was  tuned  to  attain  better  agreement  between  radar  and 
servo  winds  [Buonsanto  et  al.,  1989].  This  empirical 
tuning  was  done  with  the  help  of  radar-deduced  winds. 

There  are  two  objectives  of  the  investigation  reported 
here.  The  primary  objective  is  to  develop  a  USU  servo 
model  by  finding  the  physical  basis  of  servo  model  tuning 
and  explain/reduce  the  differences  found  in  the  FLIP  and 
servo  model  winds.  A  secondary  objective  is  to  show  how 
the  USU  servo  model  is  applied  to  different  sources  of  hmF2 
data  and  indicate  the  significance  of  this  application. 
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1.3  Objectives  of  the  Investigation 

The  main  objective  of  this  investigation  is  to  develop 
the  USU  servo  model  on  research  originally  published  by 
Rishbeth  [Rishbeth,  1967;  Rishbeth  et  al . ,  1978]  and 
[Buonsanto  et  al.,  1989],  which  calculates  mid-latitude 
meridional  neutral  winds  from  the  height  of  the  F2-peak. 

The  USU  servo  model  runs  on  a  personal  computer  (PC)  and 
requires  no  tuning.  This  enables  greater  flexibility  in 
servo  wind  calculations:  One  can  calculate  USU  servo  winds 
at  any  mid-latitude  location  under  different  ionospheric 
conditions,  without  regard  to  tuning. 

The  USU  servo  model  uses  the  Mass  Spectrometer  Inco¬ 
herent  Scatter  (MSIS-86)  model  [Hedin,  1987],  to  represent 
the  neutral  atmosphere.  The  USU  servo  model  uses  typical 
mid-latitude  estimates  of  ion  temperatures  but  is  easily 
adapted  to  use  measured  ion  temperatures.  The  USU  servo 
model  does  not  need  electron  temperature,  which  is  required 
to  calculate  the  diffusion  coefficient  in  previous  servo 
models  [Buonsanto,  1986,  1990;  Buonsanto  et  al.,  1989, 

1990;  Melendez-Alvira,  1990;  Rishbeth  and  Edwards,  1989]. 

A  second  objective  is  to  apply  the  USU  servo  model  to 
three  sources  of  hmF2  data  and  discuss  or  compare  results. 
One  source  is  the  International  Reference  Ionosphere  (IRI) 
[Bilitza  et  al.,  1987],  which  models  monthly  averages  of 
hmF2  data  for  magnetically  quiet  conditions  in  the  non- 
auroral  ionosphere.  The  main  advantage  in  using  IRI  data 
is  the  ability  to  calculate  global,  climatological  neutral 
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winds,  for  any  mid-latitude  location,  under  different  solar 
cycle  conditions. 

A  second  source  is  scaled  ionosonde  critical  frequen¬ 
cies,  which  are  used  to  determine  hmF2  values  by  an  empiri¬ 
cal  equation  developed  by  Bradley  and  Dudeney  [1973].  In 
this  study,  USU  servo  model  winds  are  calculated  from 
ionosondes  operating  during  January,  1988  (GISMOS-88  cam¬ 
paign)  .  The  ability  to  use  ionosonde  data  greatly  enhances 
global  observational  coverage.  Approximately  126  ionosonde 
stations  are  operational  and  contribute  to  the  World  Data 
Center  providing  geographically  dispersed  observational 
data  regardless  of  tropospheric  weather  conditions  [Conk- 
right  and  Brophy,  1982]. 

The  third  source  of  hmF2  data  is  IS  radar  measurements 
from  Millstone  Hill  during  the  equinox  transition  study 
(ETS)  of  September,  1984.  Two  geomagnetic  storms  occurred 
during  this  period  and  prompted  a  comparative  study  of 
FLIP,  servo  and  radar  winds  by  Buonsanto  et  al.  [1989].  By 
using  this  same  ETS  period,  the  USU  servo  model  joins  this 
comparison. 

This  investigation  reveals  the  physical  sources  of  the 
differences  between  the  Miller  et  al.  [1986]  FLIP  winds  and 
the  Buonsanto  et  al.  [1989]  servo  winds.  The  sources  stem 
from  the  theoritical  assumption  of  diffusive  equilibrium. 
This  study  also  shows  wind  results  from  different  hmF2 
sources  and  indicates  the  value  in  using  a  servo  model  for 
global  meridional  neutral  wind  calculations. 
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CHAPTER  II 

REVIEW  OF  LITERATURE 

2.1  Meridional  Neutral  Winds 

In  this  thesis,  meridional  neutral  wind  is  defined  as 
the  horizontal  motion  of  neutral  gas  particles  along  a 
geomagnetic  meridian.  The  distinction  of  neutral  wind  is 
used  at  high  altitudes  because  neutral  and  ionized  plasma 
gases  have  different  constraints  on  their  movements.  In 
the  F-region  of  the  ionosphere,  plasma  is  constrained  to 
move  along  geomagnetic  field  lines.  Plasma  motions  are 
driven  by  electromagnetic  forces,  neutral  winds,  plasma 
diffusion  and  thermal  expansion  and  contraction  [Rishbeth, 
1972;  Salah  and  Holt,  1974;  Wickwar  et  al.,  1984].  Meridi¬ 
onal  neutral  winds  cause  neutral-air  particles  to  collide 
with  F-region  ions,  forcing  them  to  move  along  inclined 
geomagnetic  field  lines.  Nighttime  equatorward  winds  raise 
the  F-layer,  which  maintains  ionization  by  decreasing 
chemical  loss,  while  daytime  poleward  winds  lower  the  F- 
layer  and  increase  chemical  loss  [Hargreaves,  1979;  Rish¬ 
beth,  1972;  Schunk,  1983].  Normally,  mid-latitude  neutral- 
wind-induced  plasma  motions  are  greater  in  magnitude  than 
electric  field  effects  [Hargreaves,  1979;  Rishbeth,  1972; 
Roble,  1983].  Methods  that  derive  neutral  winds  from  hmF2 
values  assume  that  electric  field  effects  are  negligible. 

Neutral  wind  is  mainly  created  by  an  uneven  heating  of 
the  Earth's  upper  atmosphere,  but  neutral  air  may  also  be 
set  into  motion  by  ion  collisions  created  from  electric 
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fields,  especially  important  at  high  latitudes  [Killeen  and 
Roble,  1984;  Roble,  1983].  The  driving  energy  is  provided 
by  heat  created  primarily  from  solar  ultraviolet  (UV)  and 
extreme  ultraviolet  (EUV)  radiation.  But  heat  energy  can 
be  significantly  enhanced  by  auroral  processes,  ion-neutral 
collisions  and  upward-propagating  waves  [Rishbeth,  1988; 
Richmond,  1983].  A  daytime  thermal  expansion  of  the  atmos¬ 
phere  creates  a  diurnal  bulge,  which  gives  rise  to  horizon¬ 
tal  and  vertical  pressure  gradients  accounting  for  the 
solar  driven  neutral  wind  [Rishbeth,  1972].  Although  the 
vertical  component  of  the  pressure  gradient  is  the  larger, 
it  is  almost  completely  balanced  by  gravity.  This  is  the 
foundation  of  hydrostatic  balance  and  is  an  excellent 
approximation,  especially  when  considering  large  scale 
flows  [Holton,  1979;  Richmond,  1983].  As  a  consequence, 
the  primary  driver  of  the  neutral  wind  is  in  the  horizontal 
plane. 

The  acceleration  of  horizontal  neutral  winds  is  sub¬ 
ject  to  four  types  of  momentum  [Rishbeth  and  Garriott, 

1969;  Rishbeth,  1972].  The  four  forces  are  caused  by 
pressure  gradients,  the  Earth's  rotation,  ion  drag  and 
viscous  drag.  To  best  describe  these  forces,  we  first 
introduce  the  equation  of  motion,  represented  as 

dU  n 

—  =  VP  -  2fl  x  D  -  vni<u  “  v)  +  — p"  V2U  +  g  (1) 

(a)  (b)  (c)  (d)  (e)  (f) 

where  the  terms  are  (a)  acceleration,  (b)  pressure  gradi- 


8 


ent,  (c)  Coriolis  acceleration,  (d)  ion  drag,  (e)  viscous 
drag  and  (f)  gravity.  For  the  parameters  of  Equation  (1) , 

U  is  the  neutral  wind  velocity,  P  is  pressure,  n  is  the 
Earth’s  angular  velocity,  vn^  is  the  neutral-ion  collision 
frequency,  V  is  the  ion  drift  velocity,  n  is  the  dynamic 
viscosity,  P  is  neutral  mass  density,  g  is  gravity  and  the 
ratio  n/p  is  the  kinematic  viscosity. 

As  stated,  gravity  and  the  vertical  pressure  gradient 
act  to  cancel  each  other  and  neutral  winds  tend  to  move 
horizontally.  As  a  result,  vertical  gradients  of  gas 
composition  and  temperature  are  much  greater  than  horizon¬ 
tal  gradients  and  horizontal  wind  speeds  are  usually  sever¬ 
al  orders  of  magnitude  greater  than  vertical  wind  speeds 
[Richmond,  1983;  Rishbeth,  1972].  Given  these  conditions, 
the  acceleration,  pressure  gradient,  Coriolis  force  and  ion 
drag  tend  to  act  in  the  horizontal,  while  viscous  drag  only 
appears  under  the  presence  of  vertical  wind  shears. 

As  we  are  on  the  rotating  Earth,  we  naturally  prefer 
to  use  a  reference  frame  which  is  fixed  with  respect  to  the 
Earth.  The  physical  laws  of  motion  require  an  inertial 
reference  frame,  such  that  an  object  set  in  uniform  motion 
will  remain  in  uniform  motion  in  the  absence  of  any  forces 
(Newton’s  first  law  of  motion).  Since  a  rotating  frame  is 
a  non-inertial  reference  frame,  two  "apparent"  forces  must 
be  added  to  satisfy  the  equations  of  motion.  These  appar¬ 
ent  forces  are  called  Coriolis  and  centrifugal  accelera¬ 
tions.  Centrifugal  acceleration  is  needed  to  offset  cen- 
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tripetal  acceleration,  which  acts  toward  the  axis  of  rota¬ 
tion  as  viewed  from  a  fixed  reference  frame  in  space.  The 
centrifugal  acceleration  is  an  inertial  reaction  needed  to 
balance  the  rotating  coordinate  acceleration  as  observed  on 
a  stationary  point  on  the  Earth.  The  centrifugal  force 
partly  balances  the  gravitational  force,  and  it  is  thus 
convenient  to  combine  the  effects  of  gravitational  and 
centrifugal  accelerations  into  a  term  defined  as  gravity. 

When  the  neutral  wind  appears  as  to  move  along  a 
straight  path  from  a  fixed  point  in  space  (an  inertial 
reference  frame) ,  it  is  observed  from  the  rotating  Earth  to 
follow  a  curved  path.  The  apparent  force,  which  acts  to 
deflect  the  neutral  wind  from  its  inertial  straight-line 
path,  is  called  the  Coriolis  force.  The  curved  path  is  in 
a  direction  opposite  to  the  direction  of  coordinate  rota¬ 
tion.  An  object  with  a  horizontal  relative  velocity  is 
deflected  to  the  right  in  the  northern  hemisphere  and  to 
the  left  in  the  southern  hemisphere. 

The  driving  mechanism  behind  the  pressure  gradient 
force  is  differential  heating.  Solar  forcing  is  important 
throughout  the  atmosphere,  but  at  high  latitudes  signifi¬ 
cant  heat  is  also  produced  by  auroral  particle  precipita¬ 
tion  and  ion-neutral  collisions  associated  with  the  polar 
cap  convective  cells  [Schunk,  1983].  At  mid-latitudes, 
solar  heating  is  normally  the  main  mechanism  in  pressure 
gradient  forcing.  Temperature  variations  due  to  solar 
heating  causes  the  atmosphere  to  expand  or  bulge  on  the 
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daylit  side  of  the  Earth.  This  creates  a  high-pressure 
region  referred  to  as  the  daytime  pressure  bulge  [Rishbeth, 
1972].  The  resulting  mid-latitude  daytime  poleward  and 
nighttime  equatorward  winds  have  a  significant  effect  on 
the  distribution  of  plasma  at  thermospheric  heights. 

The  effect  of  colliding  ions  on  neutral-air  particles 
is  called  ion  drag.  The  ion  drag  term  in  Equation  (1)  is 
usually  negative  at  mid-latitudes,  but  ion  drag  is  also 
capable  of  inducing  motion  to  the  neutrals  from  plasma 
drifts  created  by  electric  fields,  especially  important  in 
the  polar  regions  [Hargreaves,  1979].  At  high  latitudes, 
polar  cap  convective  cells  are  connected  to  the  magneto¬ 
sphere  and  associated  convective  electric  fields  are  en¬ 
hanced  during  ionospheric  and  geomagnetic  storms.  These 
polar  electric  fields  drive  ions,  which  collide  with  neu¬ 
tral-air  particles  having  a  significant  impact  on  the 
neutral  winds  [Schunk,  1983].  Collisions  between  electrons 
and  ions  and  between  electrons  and  neutrals  have  little 
effect  on  the  ion  drift  velocity  and  overall  neutral  air 
motion  [Kelley,  1989]. 

Physical  constraints  on  charged  ions  are  quite  differ¬ 
ent  than  on  neutral  air  movements  and  this  leads  to  some 
interesting  relations.  Neutral  air  moves  under  the  influ¬ 
ence  of  horizontal  forces,  while  at  F-region  heights,  ions 
are  constrained  to  move  along  the  geomagnetic  field  lines 
and  are  influenced  mainly  by  electric  fields  and  neutral 
air  collisions.  As  a  simple  analogy,  consider  ions  as 
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beads  "tied"  to  an  inclined  geomagnetic  string  acting  as  a 
brake  against  the  neutral  wind.  If  set  into  motion  by 
electric  fields,  the  ions  move  perpendicular  to  the  field 
lines  and  take  a  more  active  role  in  inducing  horizontal 
neutral  wind  movement. 

The  most  significant  cause  of  ion  drag  at  mid  and  low- 
latitudes  is  due  to  F-region  electric  fields  produced  in 
the  E-region  [Hargreaves,  1979;  Rishbeth,  1979].  E-region 
neutral  winds  can  blow  ions  across  the  magnetic  field  lines 
at  heights  up  to  about  150  km,  while  neutral  winds  have 
little  cross  magnetic  effect  on  electrons  above  about  70  km 
[Rishbeth,  1988].  Neutral  wind  effects  on  electron  and  ion 
motions  perpendicular  to  the  magnetic  field  lines  depend  on 
the  ratio  of  collision  frequency  to  gyrofrequency .  Gyro- 
frequency  is  inversely  proportional  to  particle  mass,  while 
collision  frequency  decreases  exponentially  with  height. 
Therefore,  neutral  winds  blow  ions  across  geomagnetic  field 
lines  at  greater  heights  than  it  does  electrons,  so  that 
ions  and  electrons  become  separated.  Bulk  neutrality  is 
maintained,  but  a  potential  difference  is  set  up  between 
ions  and  electrons,  which  is  constantly  changing  as  wind 
velocity  and  ion  concentration  change.  Neutralizing  cur¬ 
rents  are  produced  in  the  E-region.  Geomagnetic  field 
lines  act  as  equipotential  surfaces,  due  to  the  high  longi¬ 
tudinal  conductivity  along  the  field  lines  as  compared  to 
the  smaller  transverse  conductivities.  Since  the  potential 
drop  along  a  field  line  is  small,  the  electric  field  pro- 
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duced  in  the  E-region  is  reproduced  in  the  F-region.  The 
zonal  electric  field  in  the  F-region  then  acts  to  raise  or 
lower  the  F2-peak  layer  height,  where  the  ion  drift  veloci¬ 
ty  becomes  V  =  (E  x  B)/B* .  An  eastward  electric  field 
raises  hmF2  and  a  westward  field  lowers  hmF2 .  In  this  way, 
the  E-region  is  considered  an  atmospheric  dynamo  and  the  F- 
region  is  analogous  to  a  motor  [Hargreaves,  1979;  Kelley, 
1989;  Ratcliffe,  1972].  Normally,  dynamo-induced  ion 
drifts  in  the  mid-latitude  F-region  are  less  than  30  m  s'1 
[Richmond,  1983;  Rishbeth,  1972].  Miller  et  al.  [1987] 
showed  that  for  quiet  or  moderate  geomagnetic  conditions, 
electric  field  induced  drifts  are  usually  less  than  the 
overall  statistical  uncertainty  in  the  calculation  of 
neutral  winds  from  hmF2. 

The  viscosity  of  a  fluid  is  a  measure  of  its  resist¬ 
ance  to  flow.  The  greater  the  resistance,  the  greater  the 
viscosity.  Viscous  forces  appear  only  under  the  presence 
of  wind  shears.  Normally,  only  vertical  wind  shears  are 
important:  Vertical  gradients  dominate  over  horizontal 
gradients  and  horizontal  winds  are  usually  much  greater 
than  vertical  winds.  Dynamic  viscosity  depends  on  tempera¬ 
ture.  Since  F2-region  temperature  gradients  are  small,  we 
can  treat  the  dynamic  viscosity  as  a  constant  with  respect 
to  height  [Rishbeth,  1972].  Kinematic  viscosity  is  in¬ 
versely  proportional  to  neutral  mass  density,  so  that 
viscosity  is  large  and  important  at  great  heights.  The 
interchange  of  the  horizontal  momentum  in  the  vertical 
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direction  is  manifested  as  viscosity,  so  that  momentum 
tends  to  diffuse  in  such  a  way  as  co  smooth  out  the  verti¬ 
cal-height  profile.  At  about  250  km,  viscosity  acts  to 
smooth  out  velocity  profiles  in  less  than  an  hour,  while  at 
130  km,  viscous  forces  take  about  a  day  to  smooth  the 
velocity  profiles  [Richmond,  1983].  For  this  reason  we  can 
assume  that  the  neutral  wind  calculated  at  the  F2-peak  is 
constant  within  the  F-region. 

The  relative  importance  of  the  forcing  terms  in  the 
F2-region  is  variable  and  depends  on  diurnal,  geographical, 
seasonal,  solar  cycle  and  active  solar/geomagnetic  condi¬ 
tions.  Solar  activity  causes  the  greatest  variability, 
creating  geomagnetic  and  ionospheric  disturbances,  that 
affect  the  entire  thermosphere,  especially  at  high  lati¬ 
tudes.  At  solar  cycle  maximum,  increased  uv  and  EUV  heat¬ 
ing  raises  thermospheric  temperatures  by  hundreds  of  de¬ 
grees  and  changes  the  entire  dynamic  and  chemical  structure 
of  the  thermosphere  [Roble,  1983;  Rishbeth,  1986].  Al¬ 
though  the  momentum  forcing  terms  do  vary,  we  can  make  some 
general  observations  about  their  relative  significance. 

During  geomagnetically  active  periods  neutral  wind 
circulation  patterns  are  severely  disturbed  [Burnside  et 
al.,  1991].  Rapidly  changing  conditions  give  rise  to 
differential  horizontal  momentums,  so  that  viscous  forces 
maximize  to  smooth  out  resulting  vertical  wind  shears 
[Richmond,  1983].  Viscosity  then  acts  as  a  governing  force 
on  the  thermosphere,  continually  trying  to  restore  the 
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neutral  wind  back  to  quiet,  steady-state  conditions. 

At  steady  state,  we  may  neglect  acceleration  and 
viscosity  and  a  force  balance  relationship  occurs  between 
the  pressure  gradient,  ion  drag  and  Coriolis  forces.  In 
the  daytime  F2-region,  ion  drag  is  large  compared  to  Corio¬ 
lis  force,  so  that  the  neutral  wind  blows  almost  parallel 
to  the  pressure  gradient  [Rishbeth,  1988].  At  night, 
ionization  is  near  zero  and  ion  drag  is  smaller.  Nighttime 
neutral  wind  speeds  are  greater,  making  the  Coriolis  force 
more  important.  The  force  balance  relationship  at  night, 
then  tends  to  cause  neutral  air  to  flow  at  greater  angles 
to  the  pressure  gradient. 

Figure  1  shows  how  the  relationship  between  forces 
works  in  the  northern  hemisphere.  Ion  drag  acts  in  the 
opposite  direction  of  the  neutral  wind  velocity  and  Corio¬ 
lis  acts  to  the  right  of  the  wind.  The  pressure  gradient 
force  acts  from  high  to  low  pressure. 

The  force  balance  is  affected  when  pressure  gradients 
or  ion  drag  change.  One  major  effect  comes  from  solar 
cycle  variations,  which  significantly  change  pressure 
gradients  and  ion  drag.  Increases  in  the  pressure  gradient 
and  ion  drag  momentums  are  directly  related  to  increases  in 
solar  activity  [Rishbeth,  1972].  The  study  of  changes  in 
neutral  wind  force-balance  relationships  resulting  from 
changes  in  the  solar  cycle  can  lead  to  an  understanding  of 
thermospheric  circulation  patterns. 
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Day 


Night 


Fig.  1.  Day/Night  force  balance  relationship  between  the 
pressure  gradient,  ion  drag  and  Coriolis  forces  in  the 
northern  hemisphere.  U  depicts  the  horizontal  neutral 
wind.  Adapted  from  Holton  [1979]. 
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2.2  The  F2-Peak  as  a  Servo  Mechanism 

The  F2-peak  marks  the  level  of  peak  electron  density 
concentration  and  occurs  at  the  height  where  photochemical 
processes  give  way  to  gravitational  or  diffusive  control 
[Rishbeth,  1988].  Below  the  peak,  the  electron  density 
increases  with  height  because  the  loss  rate  decreases  more 
quickly  than  the  production  rate  and  neutral  densities  are 
large  enough  to  make  plasma  diffusion  of  little  importance 
as  a  loss  mechanism  [Hargreaves,  1979].  Above  the  peak, 
the  atmospheric  density  is  so  low  that  chemical  loss  gives 
way  to  loss  by  diffusive  vertical  transport.  Electron 
density  behavior  at  the  peak  is  governed  by  a  balance  of 
production,  loss  and  a  downward  diffusion  of  plasma,  where 
both  chemical  loss  and  diffusion  are  about  equally  impor¬ 
tant  [Rishbeth  and  Garriott,  1969;  Schunk,  1983]. 

Ions  and  electrons  are  held  together  by  electrostatic 
forces,  which  maintain  a  total  bulk  neutrality  [Richmond, 
1983].  In  the  F-region,  plasma  is  constrained  to  flow 
along  the  geomagnetic  field  lines,  where  meridional  neutral 
winds  or  zonal  electric  fields  can  induce  plasma  motions 
[Schunk,  1983]. 

If  the  F2-layer  is  not  disturbed  by  a  neutral  wind  or 
electric  field,  then  the  F2-peak  lies  at  a  balance  height, 
which  we  will  call  hQ.  If  the  peak  is  shifted  in  altitude, 
there  is  a  response  to  the  change  in  chemical  and  diffusive 
equilibrium  and  the  layer  seeks  to  find  a  new  equilibrium. 
In  this  way,  the  F2-peak  acts  as  a  servo  mechanism,  where 
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equilibrium  is  disturbed  by  a  vertical  drift  treated  as  a 
time-varying  perturbation  [Rishbeth,  1967;  Rishbeth  et  al., 
1978]  . 

To  derive  the  meridional  neutral  wind,  the  USU  servo 
model  first  calculates  hQ  from  a  balance  of  chemical  and 
diffusive  loss  of  plasma  at  the  F2-peak  with  no  applied 
drift.  Measured  or  empirically  modeled  F2-peak  layer 
heights  defined  as  hmF2  or  are  input  into  the  USU  servo 

model  and  applied  plasma  drifts  are  calculated  from  a 
difference  between  h^j  and  hQ.  Meridional  neutral  winds  are 
then  calculated  from  the  plasma  drift  values. 

Figure  2  shows  how  the  F2-layer  may  be  shifted  in 
altitude  by  an  applied  force.  In  Figure  2,  electron  densi¬ 
ty  profiles  are  generated  by  the  International  Reference 
Ionosphere  (IRI-86) ,  for  a  geomagnetically  quiet  period, 
with  low  solar  activity.  The  profiles  are  calculated  at 
Petropavlovsk  on  the  Kamchatka  Peninsula  in  January,  at 
1500  and  2200  local  time.  Solid  lines  show  IRI  hm  and 
dashed  lines  show  hQ  calculated  by  the  USU  servo  model  for 
each  profile.  Since  hQ  occurs  where  there  is  no  applied 
force  to  disturb  the  F2-peak,  then  any  height  difference 
must  be  attributed  to  either  a  meridional  neutral  wind, 
zonal  electric  field  or  a  combination  of  the  two.  If  we 
assume  the  electric  field  is  negligibly  small,  then  the  USU 
servo  model  takes  advantage  of  equilibrium  changes  by 
calculating  the  required  meridional  component  of  the  neu¬ 
tral  wind.  For  the  example  of  Figure  2,  calculated  values 
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ELECTRON  NUMBER  DENSITY  (PER  CM3) 


Fig.  2.  Electron  density  profiles  at  Petropavlovsk  gener¬ 
ated  with  IRI  data  for  low  solar  flux  during  winter  sol¬ 
stice.  Dashed  lines  indicate  each  profile's  balance  height 
computed  by  the  USU  servo  model.  A  computed  40  m  s'1  pole- 
ward  wind  lowers  the  daytime  F2-peak,  while  a  135  m  s'1- 
equatorward  wind  raises  the  nighttime  F2-peak. 
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of  a  daytime  poleward  wind  of  40  m  s'1  and  a  nighttime 
equatorward  wind  of  135  ms*1  are  needed  to  shift  the  F2- 
peak  layer. 

Since  we  assume  the  total  plasma  drift  causing  the  F2- 
layer  to  move  from  hQ  to  h^  is  from  the  neutral  wind,  it  is 
important  to  know  how  much  error  is  introduced  by  neglect¬ 
ing  electric  field  effects.  This  same  major  assumption  is 
used  by  any  method  deriving  neutral  winds  from  the  height 
of  the  F2-peak.  Miller  et  al.  [1987]  made  a  comparison 
study  on  the  errors  introduced  when  the  effect  of  electric 
fields  are  neglected.  By  using  the  FLIP  model,  neutral 
winds  were  derived  from  F2-peak  heights  for  quiet  and 
moderate  geomagnetic  conditions  over  Millstone  Hill  and 
Ottawa.  Coordinated  measurements  using  an  IS  radar  at 
Millstone  Hill  and  an  ionosonde  over  Ottawa,  showed  reason¬ 
able  agreement  in  both  wind  speed  and  direction,  as  well  as 
time  of  occurrence.  Electric  field  effects  produced  errors 
less  than  the  statistical  uncertainties  in  the  model  and 
measurements.  Buonsanto  et  al.  [1989]  made  a  similar  study 
using  Millstone  Hill  radar  data  and  came  up  with  the  same 
conclusions,  but  added  that  for  active  geomagnetic  periods 
electric  fields  should  not  be  neglected. 

Understanding  F2 -layer  behavior  resulting  from  iono¬ 
spheric  processes  requires  a  mathematical  basis  beginning 
with  mass  continuity  [Rishbeth,  1986] .  The  Rishbeth, 
Ganguly  and  Walker  derivation  of  the  servo  model  [Rishbeth 
et  al.,  1978],  herein  called  the  RGW  servo  model,  begins 
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with  the  F-layer  continuity  equation.  A  review  of  this 
equation  explains  the  reason  for  the  F2-peak  and  leads  to 
an  understanding  and  derivation  of  servo  relationships. 

The  F-layer  continuity  equation  for  electron/ion  density  N 
is 


dN 

at 


3  (NV7) 
3z 


(2) 


At  the  F2-peak,  where  gas  densities  are  low,  chemical 
loss  is  linear  with  N  and  is  controlled  by  charge  exchange 
reactions  of  the  major  ion  0+  with  N2  and  02  [Rishbeth, 
1986] .  The  linear  loss  coefficient  /9  is  inversely  propor¬ 
tional  to  the  N2  (or  02 )  concentration  [Hargreaves,  1979]. 
Normally  [N2 ]/[02 ]  «  10  in  the  F2-region,  but  loss  reac¬ 
tions  with  02  are  «  10  times  faster  than  with  N2  and  both 
reactions  are  considered  equally  important  [Rishbeth, 

1986].  The  F2-layer  linear  loss  coefficient  then  takes  the 
form  /3  =  kjJOj  ]  +  k2[N2  ],  where  k^.  and  k2  are  chemical  rate 
reactions  given  in  Table  1.  General  values  of  (3  are  on  the 
order  of  1  x  10'3  s"1  by  day  and  1  x  10‘4  s'1  at  night 
[Rishbeth,  1986]. 

Horizontal  motion  contributes  little  to  Equation  (2) , 
so  the  transport  term  includes  only  the  vertical  plasma 
velocity  V2  [Rishbeth,  1986].  The  value  of  Vz  includes 
both  plasma  diffusion  and  vertical  drift  from  an  applied 
wind  or  electric  field.  In  ionospheric  theory,  it  is 
sometimes  preferable  to  define  a  reduced  height  in  terms  of 
the  major  ionizable  gas,  atomic  oxygen  [Rishbeth,  1986]. 
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TABLE  1.  Ion  Chemistry  and  Reaction  Rates, 

Adapted  from  Schunk  [1983] 


Tr  =  (Tn  +  T^)/2  —  Effective  temperature,  where  Tn  and 
are  the  neutral  and  ion  temperatures  in  K. 

We  then  assume  the  scale  height  H  is  of  atomic  oxygen, 
which  is  constant  with  height,  since  the  F2-layer  is  nearly 
isothermal.  H  is  equal  to  kbTn/mg  =  RTn/Mg ,  where  kb  is 
Boltzmann's  constant,  R  is  the  universal  gas  constant,  m  is 
particle  mass,  M  is  molar  mass,  Tn  is  the  neutral  air 
temperature  and  g  is  gravity  [Hess,  1979;  Hargreaves, 

1979].  Reduced  height  in  the  F-region  is  z  =  (h  -  hr)/H 
and  is  measured  in  units  of  H,  where  h  is  the  altitude  of 
interest  and  hr  is  a  reference  height  at  a  fixed  boundary 
within  the  F2-layer.  Equilibrium  between  sources  and  sinks 
of  ionospheric  plasma  are  normally  in  balance,  except  at 
sunrise  and  sunset  [Ratcliffe,  1972].  The  servo  model  uses 
equilibrium  conditions  to  describe  F2-peak  behavior  through 
a  balance  of  terms  in  the  F-layer  continuity  equation. 
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2.3  Chemical  Loss  in  the  F2-reaion 

The  F2-layer  is  well  above  the  level  of  peak  solar 
ionization  [Hargreaves,  1979;  Rishbeth  et  al.,  1978].  This 
means  the  ion  production  rate  q  depends  on  the  photochemis¬ 
try  of  atomic  oxygen  and  is  proportional  to  [0]  and  thus 
proportional  to  exp(-z) .  As  stated,  p  depends  on  charge 
exchange  reactions  of  0+  ions  with  N2  and  02  and  is  propor¬ 
tional  to  [N2 ]  and  [02 ]  and  thus  proportional  to  exp(-kz) , 
where  k  is  the  loss  scale  height  factor.  The  use  of  the 
factor  k  allows  for  a  scale  height  ratio  adjustment  to  z, 
which  is  measured  in  units  of  atomic  oxygen  scale  height  H. 
Since  z  is  measured  in  terms  of  H,  P  can  be  evaluated  at 
any  height  in  relation  to  exp(-kz) ,  rather  than  use  differ¬ 
ently  based  scale  heights  in  N2  or  02 .  To  evaluate  k, 
consider  Equation  (2)  at  equilibrium,  so  that  the  electron 
density  Nm  a  qn t/Pm,  where  m  denotes  the  value  at  the  F2- 
peak.  The  loss  scale  height  factor  k  inversely  corresponds 
to  the  ratio  of  N2  or  02  scale  height  to  0  scale  height. 
This  becomes  a  mass  ratio  relationship  where  1.75  (28/16) 
corresponds  to  N2  loss  and  2  (32/16)  corresponds  to  02 
loss.  The  value  of  k  varies,  depending  on  loss  reactions 
of  0+  with  N2  and  02 .  The  chemical  loss  process  is  in  two 
steps,  where  the  first  step  dominates  in  the  F2 -region  and 
the  second  step  involves  dissociative  recombination  reac¬ 
tions  of  molecular  ions  with  electrons  [Hargreaves,  1979; 
Schunk,  1983;  Rees,  1989].  The  controlling  first  step 


reactions  are 
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+  02  -*■  02  +  0,  k^_ 

(3) 

+  N2  -  N0+  +  0,  k2 

(4) 

where  k-^  and  k2  are  chemical  rate  reactions  (see  Table  1)  . 
The  reaction  rate  k2  given  in  Table  1  does  not  take  into 
account  the  strong  dependence  on  the  degree  of  vibrational 
excitation  of  N2 .  Richards  and  Torr  [1986]  have  shown  that 
this  effect  is  important  at  solar  maximum  in  the  summer  and 
could  be  important  during  magnetically  active  conditions. 
When  both  N2  and  02  are  considered  equally  important,  an 
average  constant  for  k  may  be  taken  as  1.875  [Buonsanto  et 
al.,  1989].  A  better  assumption  of  the  loss  scale  height 
factor  is  to  evaluate  k  as  a  linear  interpolation  between 
loss  by  N2  and  02 .  This  improves  the  value  of  k  since 
molecular  concentrations  of  nitrogen  and  oxygen  vary  as 
well  as  reaction  rates,  which  both  depend  on  temperature 
and  density  [Schunk,  1983].  This  improved  version  of  k  is 
used  in  the  USU  servo  model  and  is  presented  in  Chapter 
III,  Section  3.1. 

?.4  Blasaa-Bifiasian 

The  concentration  of  molecular  nitrogen  and  oxygen 
decreases  with  height  more  quickly  than  atomic  oxygen,  so 
that  atomic  oxygen  dominates  the  F2-region.  Plasma  in  the 
F2-region  is  mainly  of  0+  ions  and  electrons.  Since  the 
production  of  0+  ions  depends  on  [0]  and  chemical  loss 
depends  on  [N2 ]  and  [02 ] ,  then  the  production  of  electrons 
would  increase  indefinitely  with  height  if  there  is  only 
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loss  by  chemistry.  The  F2-peak  exists  because  vertical 
transport  by  diffusion  becomes  important  as  gas  concentra¬ 
tions  decrease. 

Diffusion  of  plasma  depends  mainly  on  the  frequency  of 
collisions  between  the  ions  and  the  neutrals.  We  neglect 
electron  collisions  because  of  their  very  small  mass. 

Plasma  diffusion  is  inversely  proportional  to  the  collision 
frequencies  of  0+  moving  through  0,  N2  and  02  (see  Table  2 
for  collision  frequencies) .  Neutral  density  decreases 
exponentially  with  height,  since  the  F2-region  is  nearly 
isothermal  [Richmond,  1983].  This  means  an  exponential 
decrease  of  collisions,  which  indicates  diffusion  is  pro¬ 
portional  to  exp(z)  and  becomes  the  main  loss  mechanism 
above  the  peak.  At  the  peak,  time  scales  for  diffusion  and 
chemical  loss  are  approximately  equal.  For  an  inclined 
geomagnetic  field  at  mid-latitudes,  these  time  scales  are 
[Schunk,  1983] 


rd  = 


Jbl 


Dp0  sinJI 


Pc 


=  T, 


(5) 


where  Hp  is  the  plasma  scale  height,  I  is  the  dip  angle 
between  the  magnetic  field  and  the  horizontal  plane  and  the 
subscript  o  denotes  any  quantity  at  the  balance  height.  D 
is  the  plasma  diffusion  coefficient  (also  called  the  ambi- 
polar  diffusion  coefficient) . 

A  plasma  diffuses  along  magnetic  field  lines  in  an 
ambipolar  state  called  ambipolar  diffusion,  which  may  be 
explained  in  terms  of  gravity  and  an  electric  field 
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TABLE  2.  Collision  Frequencies  of  0+  Ions  Through  0,  N2 
and  02  Neutrals,  Adapted  from  Banks 
and  Kockarts  [1973] 


Species 


0+  -  0: 
0+  -  N2  : 
0+  -  02  : 


[n],  cm'3;  vin,  s'1 


3.67  x  10'11  [0]T J5 


6.82  x  10'10  [N2  ] 
6.64  x  10'10  [02  ] 


1  -  0.0278  In (T  ) 


vin  =  (°+  °)  +  (0+  -  N2  )  +  (0+  -»  02  ) 

Tr  =  (Tn  +  T^)/2  —  Effective  temperature,  where  Tn  and  T^ 
are  the  neutral  and  ion  temperatures  in  K.  11 

response.  In  a  plasma,  gravity  acts  to  pull  the  heavier 
ions  away  from  the  electrons.  A  restoring  electric  field 
sets  up  to  lift  the  ions  and  depress  the  electrons.  This 
causes  the  plasma  to  behave  as  thougu  each  partic1 e  has  a 
mass  of  mp  =  (m^  +  me)/2  ~  mj/2.  The  resulting  scale 
height  of  the  plasma  is  Hp  «  2kbTp/ (m^g) ,  where  the  plasma 
temperature  Tp  =  (T^  +  Te)/2.  The  ambipolar  diffusion 
coefficient  is  Dp  =  kbTp/  =  <?Hp/v^n.  At  diffusive 

equilibrium,  T^  =  Te  =  Tn  and  the  plasma  scale  height  and 
ambipolar  diffusion  coefficient  can  be  replaced  with  2H  and 
2D  [Hargreaves,  1979].  The  theoretical  RGW  servo  model 
makes  this  assumption  for  computational  convenience. 
Buonsanto  et  al.  [1989,  1990]  established  a  computational 
contradiction  by  using  this  same  assumption  for  the  plasma 
scale  height,  while  using  radar  measured  ion  and  electron 
temperatures  to  calculate  the  diffusion  coefficient. 
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Diffusion  of  plasma  along  a  magnetic  field  line  de¬ 
pends  on  density  and  temperature  gradients  as  well  as 
gravity,  all  acting  in  the  vertical  direction  [Schunk, 
1983],  Plasma  is  constrained  to  move  along  inclined  mag¬ 
netic  fields  at  F-region  heights.  Equation  (6)  shows  this 
plasma  diffusion  dependence  for  the  vertical  plasma  veloci¬ 
ty  V2,  introduced  in  Equation  (2)  and  given  as 


Vz  =  w  -  Dp  sin2 I 


N 


dN 

dh 


dT- 

dh 


H 


PJ 


(6) 


where  w  is  an  applied  vertical  drift  resulting  from  the 
meridional  neutral  wind  (neglecting  electric  field 
effects) .  At  the  F2-peak  dN/dh  =0.  If  we  assume  thermal 
equilibrium  between  plasma  and  neutral  temperatures  in  an 
isothermal  atmosphere,  we  then  have 


g  sin2 1 

vzm  =  w  "  wdm  sin‘ 1  =  w  - 

vm 


(7) 


where  the  subscript  m  refers  to  the  current  F2-peak  height 
and  WHm  is  the  plasma  diffusion  velocity  at  h^.  Here,  Dpm 
is  replaced  by  gllp/vm.  Equation  (7)  does  not  assume  diffu¬ 
sive  equilibrium  and  there  is  no  need  for  electron  tempera¬ 
ture.  The  RGW  servo  solution  uses  Dpm  in  its  Vzm  equation 
and  replaces  Hp  with  2H.  Equation  (7)  is  a  more  practical 
form  of  the  peak  vertical  plasma  velocity  and  is  used  in 
the  USU  servo  model. 

The  ion-neutral  collision  frequency  used  in  Equation 
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(7)  is  taken  from  Banks  and  Kockarts  [1973]  and  is  given  in 
Table  2.  It  has  recently  been  shown  that  the  uncertainty 
in  the  collision  cross-section  of  0+->0  warrants  an  increase 
of  1.7  [Burnside  et  al.,  1987].  This  1.7  increase  is  used 
in  the  USU  servo  model  for  applications  in  this  study. 

As  plasma  diffuses  along  an  inclined  geomagnetic  field 
line,  the  horizontal  component  of  the  plasma  motion  will 
induce  a  horizontal  neutral  wind.  The  time  it  takes  for 
the  plasma  to  induce  a  neutral  wind  depends  on  the  relative 
densities  of  the  ions  and  neutrals.  The  steady-state 
effectiveness  of  plasma  to  induce  motion  in  the  neutrals  is 
defined  as  plasma  diffusion  efficiency  and  depends  mainly 
on  the  magnetic  dip  angle.  Daytime  neutral  densities  are 
about  a  1000  times  greater  than  ion  densities  in  the  mid¬ 
latitude  F2-region  [Rishbeth,  1972].  Therefore,  downward 
diffusing  ions  take  about  1000  seconds  to  induce  neutral 
air  motion.  The  time  delay  is  regarded  as  a  time  constant 
rc j  for  ion-neutral  momentum  transfer.  Equation  (5)  gives 
this  time  constant  relationship  between  diffusion  and 
chemical  loss. 

Figure  3  shows  that  after  steady-state  conditions 
occur,  the  horizontal  component  of  the  plasma  diffusion 
velocity  equals  the  horizontal  neutral  wind,  which  is  Wd 
cot  I  [Rishbeth  and  Garriott,  1969].  Before  ion-neutral 
momentum  transfer  takes  place,  the  vertical  component  of 
the  plasma  diffusion  velocity  is  Wd  sin2 I.  As  conditions 
approach  steady  state,  the  induced  horizontal  neutral  wind 


Fig.  3.  Plasma  diffusion  along  an  inclined  geomagnetic 
field  line.  At  steady  state,  the  horizontal  component  of 
plasma  motion  induces  a  neutral  wind  which  equals  Wd  cot  I 
Adapted  from  Rishbeth  and  Garriott  [1969]. 


29 


accelerates  the  plasma  down  the  field  line  and  the  vertical 
component  of  the  diffusion  velocity  becomes  Wd.  Wd  is 
equal  to  the  vertical  velocity  as  would  occur  if  there  were 
no  magnetic  field. 

The  balance  height  in  the  servo  model  marks  the  level 
of  the  ionization  peak  if  no  ion-neutral  momentum  transfer 
has  taken  place  and  the  neutral  atmosphere  is  at  rest.  The 
balance  height  must  take  into  account  plasma  diffusion 
efficiency,  which  depends  mainly  on  dip  angle  as  shown  in 
Figure  3.  A  difference  in  plasma  diffusion  efficiency 
changes  the  balance  of  0+  ion  loss  between  diffusion  and 
chemistry.  The  USU  servo  model  calculates  balance  height 
by  taking  into  account  plasma  diffusion  efficiency  through 
a  physical  process,  not  an  empirical  one.  Buonsanto  et  al. 
[1989]  made  the  comment  that  "tuning"  the  servo  model 
physically  alters  the  relative  importance  of  loss  by  chemi¬ 
cal  recombination  and  diffusion.  Buonsanto  et  al.  [1989] 
empirically  tuned  the  servo  model  from  a  comparison  with 
radar  winds  and  increased  the  chemical  loss  importance  in 
the  servo  balance  height  relationship.  This  tuning  was 
actually  found  necessary  because  of  the  relatively  high  dip 
angle  at  Millstone  Hill  (71.5').  Tuning  the  servo  model  at 
Millstone  Hill  compensates  for  the  decreased  efficiency  of 
the  downward  diffusing  plasma  and  produces  better  agreement 
between  radar  and  servo  model  winds. 


2.5  Supporting  Ionospheric  Models 

The  USU  servo  model  depends  on  two  key  inputs  to 
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calculate  the  meridional  neutral  wind:  hmF2  as  a  function 
of  time  and  the  state  of  the  neutral  atmosphere.  The  Mass 
Spectrometer  Incoherent  Scatter  (MSIS-86)  neutral  atmos¬ 
phere  model  [Hedin,  1987]  empirically  describes  neutral 
temperatures  and  densities  and  is  used  in  the  USU  servo 
model.  The  height  of  the  F2-peak  is  either  measured  or 
empirically  modeled.  The  International  Reference  Iono¬ 
sphere  (IRI)  [Bilitza  et  al.,  1987]  provides  geomagnetical- 
ly  quiet,  monthly  averages  of  hmF2  data  for  the  non-auroral 
ionosphere.  IRI  hmF2  values  are  modeled  from  Comite  Con- 
sultatif  International  des  Radiocommunications  (CCIR)  peak 
density  maps.  CCIR  peak  densities  are  determined  from  a 
large  worldwide  data  base  of  ionosonde  critical  frequencies 
and  IS  radar  data  using  the  Bradley  and  Dudeney  [1973]  and 
Bilitza  et  al.  [1979]  empirical  methods.  In  using  IRI,  the 
USU  servo  model  calculates  magnetically  quiet,  climatologi¬ 
cal  neutral  winds  for  any  mid-latitude  location,  under 
different  solar  flux  conditions.  Both  MSIS  and  IRI  can  run 
on  IBM  compatible  personal  computers  and  therefore  compli¬ 
ment  the  USU  servo  model  quite  well. 

The  Bradley-Dudeney  model  is  used  to  determine  hmF2 
from  ionosonde  measurements  of  the  ionosphere  [Bradley  and 
Dudeney,  1973].  Figure  4  shows  the  Bradley  and  Dudeney 
model  ionosphere  consisting  of  parabolic  E  and  F  layers 
with  a  linear  region  of  electron  concentration  between 
them.  Three  input  parameters  are  required  by  the  Bradley 
and  Dudeney  model,  which  include  the  critical  frequencies 


Plasma  (rcqatncy 


Fig.  4.  Bradley-Dudeney  model  ionosphere.  The  ordinate 
and  abscissa  are  linear  in  height  and  frequency.  Reprinted 
by  permission  of  Pergamon  Press  from  Bradley  and  Dudeney 
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foE  and  foF2  and  the  ionospheric  transmission  factor 
M(3000)F2  =  MUF/foF2.  MUF  is  the  maximum  usable  frequency 
refracted  in  the  ionosphere,  which  is  received  at  a  dis¬ 
tance  of  3000  km.  These  three  parameters  are  routinely 
measured  on  a  24-hour  basis  by  more  than  126  geographically 
distributed  ionosonde  stations.  The  ionosonde  data  is 
stored  in  the  World  Data  Center  at  Boulder,  Colorado,  where 
it  is  made  available  to  the  scientific  community. 

The  calculation  of  hmF2  from  the  Bradley  and  Dudeney 
model  is  restricted  to  the  ratio  of  foF2  to  foE  greater 
than  1.7.  At  night,  when  the  E-region  has  decayed  away  or 
during  periods  of  sporadic  E,  the  values  of  foE  in  this 
investigation  are  approximated  from  the  IRI  model  densi- 
ties,  so  that  foE  =  9  x  10  (NmE)  1 ,  where  NmE  is  the  elec¬ 
tron  number  density  at  the  E-peak  layer  height  [Chen, 

1984] . 

Estimated  errors  of  hmF2  from  the  Bradley  and  Dudeney 
model  are  on  the  order  of  ±  10  km  [Bradley  and  Dudeney, 
1973;  Dudeney,  1983]  for  the  quiet  daytime,  mid-latitude 
ionosphere.  Nighttime  estimated  errors  are  expected  to  be 
larger,  due  to  the  estimated  values  of  foE.  A  recent  study 
by  Berkey  and  Stonehocker  [1989]  found  errors  of  less  than 
±  5  km  using  the  Bradley  and  Dudeney  method  in  27  hours  of 
quiet  data.  Another  recent  investigation  by  Sica  et  al. 
[1990]  concluded  that  for  quiet  periods,  estimates  of  hmF2 
from  ionosonde  data  may  have  uncertainties  of  ±  20  km. 
Errors  as  great  as  ±  20  km  could  have  a  significant  effect 
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on  neutral  wind  calculations.  More  work  is  needed  to 
accurately  determine  and  reduce  uncertainties  in  hmF2, 
which  makes  neutral  wind  calculations  more  viable. 

The  Field  Line  Interhemispheric-Plasma  (FLIP)  iono¬ 
spheric  model  [Richards  and  Torr,  1985,  1988]  is  used  to 
study  the  servo  model  balance  relationship.  The  FLIP  is  a 
one-dimensional  model,  which  uses  a  centered  dipole  approx¬ 
imation  for  the  geomagnetic  field  and  calculates  plasma 
densities  and  temperatures  along  an  entire  geomagnetic  flux 
tube  from  80  km  in  the  northern  hemisphere  through  the 
plasmasphere  to  80  km  in  the  southern  hemisphere.  The 
continuity  and  momentum  equations  are  solved  for  0+,  H+  and 
He+  ions  and  includes  ion-neutral  collisions  for  momentum 
transfer  in  the  E  and  F-regions.  Since  the  FLIP  is  a  cne- 
dimensional  model,  it  treats  neutral  density  as  a  known 
quantity  and  assumes  hydrostatic  balance,  adopting  MSIS  as 
its  neutral  atmosphere.  The  FLIP  model  requires  an  input 
of  the  meridional  neutral  wind,  which  it  normally  calcu¬ 
lates  by  the  Miller  et  al.  [1986]  method.  For  this  study, 
a  zero  wind  is  input  into  FLIP  so  that  balance  heights  (hQ) 
are  calculated.  Known  diurnal  patterns  of  FLIP  model  hQ 
values  are  used  to  study  servo  hQ  values,  to  make  the  usu 
servo  model  physically  more  realistic. 


34 


CHAPTER  III 
METHOD 


3.1  USU  Servo  Model  Development 

The  servo  model  was  originally  presented  by  Rishbeth 
[1967]  and  explicitly  derived  again  by  Rishbeth,  Ganguly 
and  Walker  [Rishbeth  et  al.,  1978].  The  development  here 
uses  the  original  RGW  equations  and  concepts  to  develop  a 
physically  improved  model  to  calculate  the  meridional 
neutral  winds.  The  RGW  servo  model  was  originally  de¬ 
veloped  for  a  theoretical  investigation  of  F2-peak  response 
to  applied  forces.  Since  the  investigation  was  theoreti¬ 
cal,  Rishbeth  et  al.  [1978]  assumed  diffusive  equilibrium 
and  constant  values  for  several  of  the  parameters  for 
computational  convenience.  In  order  to  use  the  servo  model 
in  a  practical  way,  i.e.  to  use  measured  or  empirical  data 
to  describe  the  dynamic  behavior  of  the  "real"  neutral 
atmosphere,  I  have  modified  some  of  the  RGW  equations 
presented  in  this  study. 

The  USU  servo  model  is  developed  from  an  outline  of 
the  following  principles,  which  describes  the  F2-layer  as 
it  responds  to  height  changes  in  its  peak  electron  density: 

1.  With  no  applied  vertical  drift,  the  F2-peak  lies 
at  a  balance  height  hQ,  where  the  plasma  diffusion  time 
scale  equals  the  chemical  loss  time  scale  of  0+  ions. 

2.  An  applied  meridional  neutral  wind  (electric  field 
effects  are  neglected)  shifts  hQ  to  a  new  level. 
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3.  At  any  instant,  the  actual  F2-peak  tries  to  main¬ 
tain  a  balance  of  diffusion  and  chemical  loss,  with  an 
exponential  time  constant. 

4.  The  loss  scale  height  factor  is  a  function  of 
time,  since  reaction  rates  and  neutral  densities  are  a 
function  of  time,  i.e.  k  =  k(t) . 

5.  The  0+  scale  height  is  a  function  of  time  and,  to 
a  much  lesser  extent,  location.  Temperature  varies  with 
time  and  gravity  varies  with  latitude  (gravity  is  computed 
at  an  F2-layer  reference  height),  i.e.  H  =  H(t,g) . 

The  USU  servo  model  interpolates  the  loss  scale  height 
factor  k  from  the  loss  of  0+  with  N2  and  02  between  chemi¬ 
cal  reaction  Equations  (3)  and  (4) .  This  linear  interpola¬ 
tion  equation  is  given  as  (see  Table  1) 


_ _ 

k^/kj,  +  [N2  ]/[02  ] 


x  0.25  +  1.75 


Figure  5  shows  mid-latitude  variations  of  calculated  k 
values  using  Equation  (8)  for  low,  medium  and  high  10.7  cm 
solar  flux.  Neutral  densities  and  reactions  are  based  on 
MSIS-86  [Hedin,  1987].  The  value  of  k  changes  mainly  on  a 
diurnal  cycle,  but  Figure  5  shows  that  k  also  depends  on 
changes  in  the  neutral  atmosphere  during  a  solar  cycle. 
Variations  in  k  would  approximate  the  same  changes  as  shown 
in  Figure  5,  which  would  reflect  changes  in  neutral  densi¬ 
ties  and  rate  reactions  as  other  ionospheric  parameters  are 
allowed  to  vary.  The  variation  of  k  of  about  1.3%,  shown 


Scale  Height  Factor, 


k/Solar  Flux  Comparison 


Fig.  5.  January  mid-latitude  variations  of  loss  scale 
height  factor  k  for  low.  medium,  and  high  solar  flux  of  70, 
130  and  190  s.f.u.  (10'22  W  m*2  Hz'1). 
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in  Figure  5,  would  result  in  a  variation  of  about  5  to  10% 
in  the  derived  neutral  winds.  Values  of  k  approximate 
1.875  [Buonsanto  et  al.,  1989],  but  Figure  5  indicates  that 
k  is  variable  and  should  be  adjusted  for  different  diurnal 
and  ionospheric  conditions. 


3.1.1  Model  Equations.  The  USU  servo  model  equations 
are  based  on  the  RGW  servo  model  [Rishbeth  et  al.,  1978]. 
Some  of  the  RGW  servo  concepts  and  solutions  are  reproduced 
here  for  clarity  and  for  further  development  and  discus¬ 
sion.  A  complete  derivation  of  the  servo  equations  are 
presented  by  Rishbeth  et  al.  [1978]  and  compliment  the  USU 
servo  model  development  presented  here. 

We  begin  the  derivation  by  multiplying  the  F-layer 
continuity  Equation  (2)  by  H  and  integrating  from  the  F2- 
peak  to  the  top  of  the  ionosphere.  For  this  evaluation, 
let  the  reduced  height  be  z  =  (h  -  hQ)/H,  where  hQ  is  the 
balance  height.  At  the  F2-peak,  z  =  zffl.  This  implies  that 
zm  =  (^m  “  hQ)/H,  so  zm  -  0  then  the  F2-peak  lies  at  the 
balance  height:  Thus  we  do  not  apply  a  drift.  The  result¬ 
ing  equation  is  then 


<00 

N  dz 
zm 


poo 

H  (q  -  0N)dz  -  *«  +  Nm(Vzm  -  H  dzm/dt)  (9) 
J  zm 


where  Nm(Vzln  -  H  dz^/dt)  is  the  relative  vertical  plasma 
flux  at  the  F2-peak  and  is  the  vertically  upward  plasma 
flux  at  the  top  of  the  ionosphere  and  equal  to  the  limiting 
value  of  NVZ.  Consider  the  following  RGW  assumptions  to 
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evaluate  Equation  (9)  : 

1.  The  electron  density  profile  above  the  F2-peak 
decays  at  the  same  rate  as  the  peak  electron  density  NmF2 , 
i.e.  the  height-integrated  ion  content  is  aHNm.  The  top¬ 
side  shape  factor  a  is  defined  as  [Melendez-Alvira,  1990] 

1  f00 

a  *  -  N ( z )  dz  (10) 

N 


The  topside  layer  remains  at  a  constant  shape  given  by  the 
Chapman  alpha  layer  as  a  =  2.82  and  is  used  in  the  USU 
servo  model. 

2.  The  F2-peak  is  well  above  the  level  of  peak  pro¬ 
duction,  so  q(z)  =  qmexp(-z). 

3.  The  linear  loss  coefficient,  0  is  proportional  to 
exp(-kz) ,  where  k  is  the  loss  scale  height  factor  depended 
on  the  dominant  ion  0+  chemical  reactions  with  N2  and  02  , 
given  by  Equation  (8) .  The  chemical  loss  rate  of  the 
electron  density  profile  above  the  peak  is  /3N(z)  « 
0mNmexP("k2> • 

We  now  integrate  Equation  (9)  and  get 


aH 


dN 


dt 


-  %H  “ 


W.H 


-  *oo  +  Nffl(V 


zm 


-  H  dzm/dt) 


(ID 


The  peak  vertical  plasma  velocity,  Vzm  in  Equation 
(11)  is  replaced  by  w  -  (g/vm)sin2I  from  Equation  (7).  By 
replacing  Vzm  with  Equation  (7)  there  is  no  need  for  elec¬ 
tron  temperature  (see  Table  2  for  collision  frequency) . 
Equation  (7)  does  not  assume  diffusive  equilibrium  and  is 
more  practical  to  use  than  the  original  RGW  servo  solution, 
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which  is  V2m  =  w  -  (Dpm/2H) sin2 I .  We  therefore  substitute 
Equation  (7)  into  (11),  divide  through  by  NmH  and  then 
rearrange  to  get 


=  qm  -  »«/H 
dt  *m 


a  dN 


N, 


m 


dt 


g  sin2 1 

v_H 

m 


w 

H 


(12) 


An  approximate  F-layer  continuity  equation  for  the 
peak  electron  density,  introduced  by  Rishbeth  et  al. 

[1978],  is  quite  useful  in  getting  a  reasonable  solution  to 
Equation  (12).  We  assume  that  any  flux  contributions  from 
the  top  of  the  ionosphere  is  evenly  distributed  over  the 
equivalent  thickness  otH  of  the  topside  layer.  This  conti¬ 
nuity  approximation  is 


*  qm  “  cVm 


Is 

aH 


(13) 


The  coefficient  c  in  Equation  (13)  is  empirical  of  order 
unity  and  is  introduced  to  take  account  of  transport  proc¬ 
esses  at  the  peak  [Rishbeth  et  al.,  1978].  This  coeffi¬ 
cient  was  first  presented  and  used  in  Rishbeth' s  servo 
model  equations  in  1967.  it  is  also  part  of  the  propor¬ 
tionality  constant  used  to  make  the  diffusion  and  chemical 
loss  time  scales  equal  at  the  balance  height. 

To  continue  the  model  development,  we  substitute 
Equation  (13)  for  dNm/dt  in  Equation  (12)  giving 


(1  -  a) 


USL 


N 


m 


kac  -  1 


0 


m 


g  sin2 1 


vmH 


w 


H 


(14) 
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Equation  (14)  is  actually  a  generalized  form  of  the  final 
servo  equation.  To  be  useful,  we  need  to  make  a  diurnal 
evaluation  of  Equation  (14). 


3.1.2  Daytime.  The  daytime  balance  relationship 
between  diffusion  and  chemistry  is  found  from  Equation 
(14),  where  w  =  0  and  dzm/dt  =  0.  Assuming  the  change  in 
peak  electron  density  is  slow  during  the  day  (dNm/dt  =  0  at 
equilibrium) ,  an  equilibrium  approximation  for  the  peak 
production  from  Equation  (13)  is 

%n  “  <=dW  +  (15) 


The  subscript  d  is  meant  to  indicate  that  cd  is  a  daytime 
derivative  of  the  diurnal  coefficient  c.  Now,  assume  ion 
production  and  chemistry  dominate  the  daytime  vertical 
plasma  flux,  so  a  balance  relationship  between  PQ  and  vQ 
from  Equation  (14)  gives 


k  g  sin*  I 
kcd_1  voH 


(16) 


When  using  measured  hmF2  data  from  ionosondes,  the 
resolution  is  too  poor  to  determine  an  adequate  layer 
height  response  with  respect  to  time.  This  implies  a 
steady  state  plasma  drift  assumption,  where  dzm/dt  =  0. 

This  assumption  does  not  greatly  affect  the  overall  calcu¬ 
lation  of  the  plasma  drift  [Buonsanto  et  al.,  1989].  In 
fact,  studies  show  w  as  mainly  linear  with  respect  to 
changes  in  layer  height,  especially  if  the  drifts  are  small 
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[Rishbeth  et  al.,  1978;  Miller  et  al.,  1986;  Buonsanto  et 
al.,  1989].  If  we  do  not  neglect  in  Equation  (14),  we 
can  substitute  Equation  (15)  for  qm  in  Equation  (14)  and 
solve  for  w  to  get  the  daytime  plasma  drift  as 


g  sin2  I 

w  =  -  +  H 

vm 


flmd  -  cd*) 
k 


+ 


-  1) 


aHN. 


m 


(17) 


We  find  values  for  and  vm  by  assuming  =  )3rexp  (-kzmr) 
and  vm  =  vrexp(-zmr),  where  zmr  =  (1^  -  hr)/H  and  hr  is  a 
convenient  reference  height  within  the  F2-region.  To  use 
Equation  (17),  we  can  neglect  the  flux  term,  but  we  need  to 
know  cd  from  the  balance  relationship  in  Equation  (16) . 

Consider  the  plasma  exchange  boundary  flux  through  the 
top  of  the  ionosphere  and  recall  that  fcgo/aH  assumes  the 
flux  is  equally  distributed  throughout  a  layer  thickness  of 
a H.  This  indicates  that  there  is  no  localized  flux  differ¬ 
ences  in  the  ionosphere/plasmasphere  boundary  allowed  in 
this  model.  This  means  w  should  be  insensitive  to 
unless  the  ratio  of  the  exchange  flux  to  Nm  is  very  large. 
This  argument  is  understood  by  examining  the  flux  term  in 
Equation  (17) .  Compare  a  typical  flux  value  used  by  Rish¬ 
beth  et  al.  [1978]  of  ^  a  1.4  x  10^  m'2  s'1  to  a  typical 
daytime  value  of  Nm  a  1  x  1012  m'3  [Hargreaves,  1979].  In 
this  case,  less  than  a  meter  per  second  is  contributed  to  w 
from  flux  effects,  leading  to  the  conclusion  that  daytime 
#3,  can  be  neglected. 

One  final  note  is  the  use  of  a  servo  analysis  in 
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finding  a  value  for  from  a  neutral  wind.  Due  to  the 
uncertainties  in  the  servo  equation  parameters  and  the 
large  daytime  insensitivity  of  w  to  a  change  in  it 

seems  doubtful  that  any  accurate  value  of  the  topside  flux 
can  be  ascertained.  Yet  an  upward/downward  flux  sense  may 
be  found.  This  type  of  study  was  successfully  done  using  a 
similar  servo  equation  by  Buonsanto  [1986],  relating  day¬ 
time  boundary  flux  comparisons  to  seasonal  change. 

If  we  solve  Equation  (16)  for  cd  and  substitute  the 
result  into  Equation  (17)  with  =  o,  our  result  becomes 
Equation  (18).  Since  we  used  Equation  (7)  for  the  replace¬ 
ment  of  Vzm  in  Equation  (11),  Equation  (18)  does  not  assume 
diffusive  equilibrium  and  does  not  need  electron  tempera¬ 
ture  for  the  calculation  of  w.  Equation  (18)  is  a  similar 
form  of  the  RGW  servo  equation  used  by  Buonsanto  et  al. 
[1989].  Equation  (19)  is  a  linear  form  of  Equation  (18), 
which  gives  results  that  are  approximately  the  same  as  the 
results  from  the  nonlinear  Equation  (18)  [Rishbeth  et  al., 
1978;  Buonsanto  et  al.,  1989].  These  equations  are 

g  sin2  I  r  -I 

w  =  -  [exp(zm)  -  exp(-kzm)J  (18) 

vo 

g  sin2  I 

w  ~  =  -  z  (k  +  1)  (19) 

vo 

where  again  zm  =  (h^  -  hQ)/H.  To  calculate  w  from  Equation 
(18)  or  (19),  we  must  find  the  balance  height  instead  of 
cd,  as  in  Equation  (17) .  From  the  balance  height,  we  can 
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calculate  vQ  and  zm  to  compute  the  plasma  drift. 

Equation  (16)  provides  the  means  to  derive  the  USU 
servo  balance  height.  It  is  convenient  to  use  an  F2- 
region  reference  height,  where  /3Q  =  0rexp(-kzor)  and  vQ  = 
vrexp(-zor)  and  where  zQr  =  (hQ  -  hr)/H.  Using  Equation 
(16)  to  relate  vQ  to  (3Q  gives  the  servo  balance  height 


^o  ^r  + 


(k  +  1) 


H/?rvr(kcy  -  1) 
kg  sin2  (I) 


where  cx  is  used  instead  of  cd  to  indicate  that  this  coef¬ 
ficient  is  variable.  The  variation  of  cx  is  mainly 
diurnal,  but  it  also  changes  with  location  or  ionospheric 
conditions.  To  calculate  cx,  we  can  use  either  of  the  two 
following  equivalent  equations 


(21) 

(22) 


We  can  use  FLIP  model  balance  heights  to  study  the  varia¬ 
tion  of  cx  and  calculate  balance  heights  for  different 
ionospheric  conditions  and  locations.  Adjusting  the  value 
of  cx  is  what  Buonsanto  et  al.  [1989,  1990]  describes  as 
tuning  the  servo  model. 

3.1.3  Nighttime.  The  nighttime  balance  relationship 
between  /3Q  and  vQ  is  different  from  daytime,  since  qm  =  0 
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in  Equation  (14).  The  same  procedures  used  to  find  the 
daytime  balance  Equation  (16)  are  used  for  a  nighttime 
balance  relationship,  with  the  result  that 


k  g  sin2  I 

kacn-1  voH 


(23) 


The  added  feature  of  the  topside  shape  factor  a  appears  at 
night  and  causes  the  night  balance  heights  to  rise  about 
half  a  scale  height  above  the  day  balance  heights.  Here, 
the  subscript  n  in  cn  denotes  night,  where  cn  is  slightly 
larger  than  c^.  The  value  of  cx  in  Equation  (20)  contains 
a  when  night  occurs,  so  that  cx  =  cd  during  the  day  and  cx 
=  acn  at  night.  Over  a  24  hour  period  a  typical  range  for 
cx  is  about  1.3  for  day  to  about  5.0  for  night. 

To  derive  a  nighttime  plasma  drift,  we  use  the  same 
procedures  as  with  the  daytime  Equation  (17) .  We  use 
nighttime  values  in  Equation  (14),  namely  qm  =  0,  to  get 


g  sin2  I  H/3_(l  -  kac„) 

w  =  — -  +  — ^  (24) 

vm  k 

If  we  solve  Equation  (23)  for  acn  and  substitute  the  result 
into  Equation  (24) ,  our  result  becomes  Equation  (18) .  This 
makes  Equations  (18)  or  (19)  the  more  convenient  forms, 
since  we  can  use  them  for  day  or  night  plasma  drift  calcu¬ 
lations.  The  servo  equations  so  far  developed  are  specifi¬ 
cally  for  either  day  or  night  and  therefore  a  method  is 
required  to  determine  the  diurnal  behavior  of  cx.  We  can 


then  calculate  the  plasma  drift  for  day  and  night  as  well 


as  the  morning  and  evening  transitions. 


3.1.4  Twilight  Assumptions.  To  find  values  of  cx  for 
the  morning  and  evening,  one  must  define  when  the  morning 
and  evening  periods  begin  and  end  in  terms  of  solar  zenith 
angle  X.  Rishbeth  [1967]  arbitrarily  chose  limits  of  the 
transitions  to  be  96°  and  82°  for  morning  and  evening. 

When  calculating  cx  from  FLIP  model  balance  heights,  it  is 
found  that  the  boundary  limits  for  X  are  closer  to  102°  and 
82°.  Another  method  in  finding  the  approximate  F-region 
transition  limits  is  to  use  the  aid  of  the  IRI  model.  By 
noting  the  time  of  decay  and  building  of  the  E-layer,  it  is 
found  that  the  limits  are  also  closer  to  102°  and  82°.  One 
other  method  will  get  an  approximate  value  for  x,  in 
relation  to  when  the  F2-region  may  "turn  on  and  off." 

Assume  the  Sun's  rays  are  straight  and  parallel,  so  a 
relationship  involving  the  Earth's  radius,  height  above  the 
Earth's  surface  and  X  is  R/(hr  +  R)  =  sin(180°  -  x) ,  where  R 
is  the  Earth's  radius  and  hr  is  a  reference  height.  From  a 
typical  F2-layer  height  of  300  km  and  R  =  6370  km,  the 
Sun's  rays  first  affect  the  F2-region  when  X  «  107°.  This 
simple  analysis  does  not  take  into  account  ionizing  effi¬ 
ciency,  absorption  cross-sections,  solar  intensity  or  any 
other  feature  which  may  affect  when  morning  begins  or 
evening  ends.  This  analysis  does  give  an  indication  of 
what  the  upper  boundary  should  be  and  tends  to  support  an 
upper  limit  of  102°  versus  96°. 
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Given  the  above  arguments  a  diurnal  pattern  is  broken 
into  four  periods: 


1. 

Morning 

102°  >  x  >  82° , 

x  decreasing; 

2. 

Day 

X  <  82  0  ; 

3. 

Evening 

82 “  <  X  <  102 

°,  X  increasing; 

4. 

Night 

X  >  102 '  and 

f 

H 

o 

• 

Mornings  and  evenings  tend  to  last  about  2  hours  at  mid¬ 
latitudes.  These  twilight  period  limits  make  it  possible 
to  describe  day  and  night  balance  height  relationships. 

The  next  problem  is  to  determine  a  discontinuity 
function  for  cx  based  on  the  defined  solar  zenith  angle 
transition  limits.  Rather  than  use  a  Chapman  function  for 
the  transitions  [Buonsanto  et  al.,  1990],  the  FLIP  model  is 
used  to  analyze  a  twilight  discontinuity  pattern  for  dif¬ 
ferent  ionospheric  conditions  and  locations.  Figure  6 
shows  a  calculated  diurnal  pattern  of  cx  values  from  FLIP 
and  USU  servo  model  balance  heights.  Transition  assump¬ 
tions,  as  given  above,  are  used  in  the  servo  calculation 
and  show  agreement  with  the  FLIP  model. 

Figure  6  also  indicates  that  morning  and  evening 
transitions  are  approximately  linear  with  time  between  day 
and  night.  The  cx  diurnal  pattern  for  morning  transitions 
are  more  abrupt  than  evening  transitions.  A  sharp  drop  in 
cx  occurs  during  the  first  hour  of  the  morning,  nearly 
reaching  a  daytime  value.  This  seems  physically  plausible 
considering  the  nighttime  decay  of  ionization.  With  the 
advent  of  morning,  production  dominates  and  the  ionosphere 


is  literally  "turned  on."  At  evening  there  is  a  more 
gradual  decay  of  the  ionosphere,  which  continues  at  night 
in  a  shape  preserving  manner  [Rishbeth  and  Garriott,  1969]. 

A  final  note  about  the  transition  periods  is  that 
equilibrium  no  longer  applies,  which  is  a  major  assumption 
in  the  derivation  of  the  servo  model.  Uncertainties  are 
largest  during  twilight  periods  due  to  rapidly  changing 
conditions.  A  20%  uncertainty  in  cx,  which  is  similar  to 
differences  shown  in  Figure  6,  results  in  about  a  12% 
uncertainty  in  the  derived  neutral  wind  speeds. 

3.1.5  Balance  Height  Adjustment.  The  empirical  value 
for  the  diurnal  coefficient  c  in  Equation  (13)  was  adjusted 
by  Buonsanto  et  al.  [1989;  1990]  to  "tune"  the  servo  model 
winds  to  Millstone  IS  radar  winds.  For  equal  loss  of  0+ 
with  02  and  N2  ,  c  is  usually  taken  from  Rishbeth  as  cd  = 
1.33  and  cn  =  1.73  [Rishbeth,  1967;  Buonsanto  et  al.,  1989, 
1990].  Rishbeth  [1967]  used  empirically  derived  values  of 
c  for  a  dip  angle  of  45°,  which  were  found  to  be  appropri¬ 
ate  for  Arecibo  (I  =  48.5°)  [Rishbeth,  1967;  Rishbeth  et 
al.,  1978;  Ganguly  et  al.,  1980;  Melendez-Alvira,  1990]. 

Tuning  the  servo  model  by  adjusting  c  physically 
alters  the  balance  height,  where  0+  loss  by  diffusion  and 
chemistry  is  equally  important.  Buonsanto  et  al.  [1989, 
1990]  found  it  necessary  to  decrease  c,  which  decreases  the 
balance  height,  making  diffusive  loss  of  0+  less  important 
at  Millstone  Hill  (I  =  71.5°).  The  reason  for  servo  model 
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tuning  stems  from  a  difference  in  plasma  diffusion  effi¬ 
ciency,  which  is  mainly  dependent  on  dip  angle  (see  Figure 
3) .  Horizontal  neutral  wind  motion  induced  by  ion-neutral 
collisions  is  depended  on  plasma  diffusion  efficiency. 
Adjusting  the  value  of  c  alters  the  diffusion  transport 
effects  in  Equation  (13),  which  leads  to  a  change  in 
balance  height.  Tuning  is  therefore  a  way  to  correct  or 
adjust  the  balance  height  for  differences  in  plasma  diffu¬ 
sion  efficiency. 

Instead  of  tuning  the  value  of  c,  the  USU  servo  model 
adjusts  hQ  based  on  the  boundary  conditions  that  during  the 
day  c  =  cd  =  1  and  at  night  c  =  cn  =  1.3,  at  a  dip  angle  of 
90*.  These  boundary  conditions  are  determined  by  using  the 
continuity  Equation  (13)  for  equilibrium  conditions.  At  I 
=  90”,  where  diffusive  ion-momentum  transfer  to  the  neu¬ 
trals  is  zero,  qm  =  0mNm  during  the  day  and  therefore  c  = 
cd  =  1.  We  assume  a  nighttime  boundary  by  a  comparative 
linear  regression  using  the  Rishbeth  [1967]  values,  where 
1.33 (.75)  =  1  and  1.73 (.75)  =  1.3.  The  value  of  hQ  is  then 
found  by  the  following  equation  [Rishbeth  et  al.,  1978] 

H  ln(sin2 I) 
ho  =  ^obc  "  1 

where  h0bC  is  the  boundary  condition  balance  height  calcu¬ 
lated  from  Equation  (20) .  In  Equation  (20) ,  the  boundary 
conditions  for  day/night  are  thus:  cx  =  cd  =  1  and  cx  =  acn 


(2.82) (1.3) . 
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Equation  (25)  is  developed  from  Equation  (18) .  In 
Equation  (18) ,  w  corresponds  to  a  steady-state  difference 
of  plasma  diffusion  velocity,  where  Wd  sin2 I  occurs  before 
ion-neutral  momentum  transfer  takes  place  and  Wd  occurs 
after  steady-state  conditions  apply,  with  an  induced  hori¬ 
zontal  neutral  wind  of  Wd  cot  I  (see  Figure  3) .  At  the 
boundary  condition  of  I  =  90s,  hQ  =  hQt)C.  The  term 
(H) ln(sin2 I)/k+l  is  negative  and  allows  for  an  increase  in 
balance  height  to  compensate  for  increased  plasma  diffusion 
efficiency  at  smaller  dip  angles.  When  we  adjust  hQ  from 
Equation  (25)  for  a  decrease  in  dip  angle  and  input  the 
resulting  hQ  back  into  Equation  (21)  or  (22)  to  calculate 
cx,  we  find  that  cx  becomes  larger  than  the  value  used  to 
calculate  hobc.  Larger  cx  values  produce  higher  hD  alti¬ 
tudes  because  diffusion  becomes  more  important. 

The  calculation  of  new  cx  values  are  not  implemented 
in  the  USU  servo  model  unless  further  adjustment  is  re¬ 
quired  to  cx,  as  will  be  discussed  in  Section  3.1.7.  If  cx 
values  are  adjusted,  then  we  calculate  the  final  hQ  values 
with  Equation  (20) . 

It  is  useful  to  compare  calculated  cx  values  with 
radar-tuned  values  or  FLIP  cx  values  for  the  same  location 
and  ionospheric  conditions.  Figure  6  is  one  example,  which 
illustrates  a  comparison  between  adjusted  servo  cx  values 
and  calculated  FLIP  cx  values. 

The  benifit  of  tuning  or  adjusting  the  servo  model 
becomes  significant  with  changes  in  dip  angle  or  large 
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changes  in  ionospheric  conditions.  Buonsanto  et  al.  [1989] 
showed  that  tuning  the  servo  model  at  Millstone  Hill  re¬ 
sulted  in  a  change  in  neutral  wind  by  as  much  as  29  m  s'1, 
which  was  about  a  35%  change  in  the  derived  wind.  By  using 
Equation  (25) ,  the  USU  servo  model  is  self  adjusting  and 
can  be  used  for  any  mid-latitude  location,  under  different 
ionospheric  conditions,  without  the  burden  of  tuning. 

3.1.6  Topside  Shape  Factor.  The  topside  ionosphere 
extends  from  the  F2-peak  to  the  ionosphere/plasmasphere 
boundary  (800  -  1000  km)  [Schunk,  1983;  Moyer,  1976].  We 
assume  the  nighttime  decay  of  the  topside  electron  density 
profile  occurs  at  the  same  rate  of  decay  as  NmF2 ,  where  a 
"shape  preserving"  solution  of  the  diffusion  equation  is 
given  as  [Rishbeth  and  Garriott,  1969] 


N(z,t)  =  Nr(z)exp(-0't) 


(26) 


We  evaluate  the  effective  chemical  decay  time  coefficient 
13'  by  using  Equation  (13)  with  qm  =  0  and  assuming  that 
chemical  loss  at  the  peak  dominates  vertical  plasma  flux 
effects,  so  that 


-cn^mNm 


Nm^J  =  Nm  ( 0 )  exP  ( -cn^mt ) 


(27) 


For  an  equivalent  rate  of  decay  for  the  topside  ionosphere 
and  NmF2 ,  it  is  apparent  from  Equations  (26)  and  (27)  that 

P'  =  cnV 

Transport  effects  from  c  in  Equation  (13)  are  related 
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to  the  topside  decay  rate.  We  can  compare  the  numerical 
values  of  cn  to  cd  by  qualitatively  comparing  nighttime  and 
daytime  decay  rates,  used  to  adjust  /3m  to  acquire  an  effec¬ 
tive  decay  coefficient  f3' .  The  effective  decay  decreases 
during  the  day,  approaching  /?m,  because  the  ion  production 
is  most  rapid  below  the  F2-peak  and  ceases  at  night.  The 
value  of  /3 '  is  increased  at  night  to  compensate  for  the 
decrease  of  (3m  at  a  higher  nighttime  balance  height,  with 
qm  =  0:  Hence  cd  is  numerically  less  than  cn. 

The  USU  servo  model  assumes  a  topside  shape  factor  a 
of  2.82,  from  the  Chapman  alpha  layer  solution  [Rishbeth  et 
al.,  1978].  Comparisons  of  cx  values  from  the  FLIP  model 
shows  that  an  assumption  of  a  =  2.82  is  reasonable.  The 
combination  of  acn  is  included  in  the  cx  calculations  from 
the  FLIP  model.  Keeping  a  apart  from  cn  is  useful,  since 
any  variability  of  cx  from  FLIP  calculations  is  attributed 
to  either  a  or  cn  for  future  study. 

The  assumption  that  vertical  plasma  flux  effects  are 
negligible  in  the  nighttime  decay  of  the  topside  ionosphere 
may  warrant  further  scrutiny.  A  recent  paper  by  Sica  and 
Schunk  [1990]  shows  that  the  ionosphere/plasmasphere  verti¬ 
cal  boundary  flux  is  important  in  finding  the  nighttime 
height  of  the  F2-peak.  If  can  appreciably  change  the 
topside  shape,  then  a  relationship  may  exist  between  ^  and 
a.  As  an  example,  after  a  geomagnetic  storm,  upward  plasma 
flux  tends  to  refill  the  depleted  plasmasphere  from  the 
topside  ionosphere  [Hargreaves,  1979;  Sica  and  Schunk, 
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1990].  Slab  thickness,  which  is  a  ratio  measurement  of  ion 
content  to  NmF2 ,  significantly  increases  for  a  few  days 
after  a  storm  and  is  an  indication  of  upward  plasma  flow 
[Hargreaves,  1979],  Slab  thickness  is  related  to  the 
topside  shape  factor,  as  indicated  in  Equation  (10).  An 
expected  increase  should  then  occur  in  a  after  a  storm 
event.  If  a  is  allowed  to  increase,  then  hQ  from  Equation 
(20)  will  increase,  implying  hmF2  should  also  increase  at 
night.  This  is  the  result  shown  by  Sica  and  Schunk  [1990], 
where  a  constant  upward  plasma  flux  raises  the  nighttime 
F2-peak  by  as  much  as  30  km,  if  holding  NmF2  constant. 

The  implication  is  that  perhaps  there  is  a  relation¬ 
ship  between  the  topside  shape,  which  can  be  measured  by  IS 
radar  and  vertical  plasma  flux,  which  is  usually  measured 
by  satellite  or  calculated  by  a  comprehensive  plasma  model. 
Still,  hQ  is  not  very  sensitive  to  an  increase  in  topside 
shape  factor,  where  a  factor  of  3  increase  in  a  is  needed 
to  increase  hQ  by  about  30  km.  The  reason  for  this  rela¬ 
tively  large  increase  in  a  is  probably  from  the  assumption 
that  NmF2  decays  at  the  same  rate  as  the  topside  electron 
density  profile,  which  is  physically  more  realistic  than 
holding  NmF2  constant.  Future  studies  with  cx  may  show  how 
nighttime  plasma  flux  changes  the  decay  rate  of  the  F2- 
peak.  This  may  then  lead  to  a  way  to  approximate  vertical 
plasma  exchange  flux  from  calculations  of  cx  =  acn  using 
ground  based  measurements. 
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3.1.7  Daytime  Steady  Rate  Assumption.  During  the  day, 
the  flux  term  in  Equation  (13)  is  normally  neglected,  since 
production  and  chemistry  dominate.  Under  steady-state 
conditions,  qm  =  cd/9mNm •  As  shown  in  Figure  6  and  again  in 
Figure  7,  the  FLIP  model  cx  value  has  a  fairly  linear 
daytime  slope.  Earlier  servo  models  assume  daytime  steady- 
state  conditions  of  cx  =  cd  =  constant.  Actually,  the 
condition  qm  =  cd^mNm  occurs  close  to  solar  noon,  where 
dNm/dt  =  0.  Before  solar  noon,  ion  production  is  greater 
than  loss  and  dNm/dt  >  0.  After  solar  noon  dNm/dt  <  0. 

Thus  >  cd0mNm  before  noon  and  qm  <  after  noon. 

To  correct  for  the  steady-state  assumption,  the  USU 
servo  model  has  an  option  of  assuming  a  daytime  steady- 
rate,  where  cd  *  cd(t)  .  The  FLIP  model  calculation  of  cx 
indicates  that  cd(t)  changes  linearly  with  time,  so  we  may 
assume  a  steady  rate.  This  means  that  qm  =  c^tJ/J^N^, 
where  dNm/dt  =  0  for  the  entire  day.  Figure  7a  shows  how  a 
steady-state  assumption  compares  with  the  FLIP  results 
during  January  at  a  mid-latitude  location.  Compare  the 
daytime  and  nighttime  deviations  between  FLIP  and  USU  servo 
cx  values  and  note  the  larger  deviation  in  day  and  night 
balance  heights.  This  is  because  daytime  cx  values  are 
close  to  unity  and  therefore  the  percent  change  in  cx 
during  the  day  are  relatively  large  compared  with  the 
percent  differences  of  the  nighttime  cx  values,  where  cx  = 
acn.  This  justifies  the  need  to  model  the  daytime  cx 
pattern  and  assume  a  steady  rate  of  change  instead  of  a 
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constant  value  as  in  a  steady-state  assumption.  Figure  7 
indicates  that  a  30%  difference  between  FLIP  and  USU  servo 
daytime  cx  vaxues  will  result  in  about  a  15  km  difference 
in  balance  height,  which  translates  to  about  a  14  m  s'1 
difference  in  neutral  wind.  To  decrease  the  FLIP/USU  servo 
difference  in  cx,  Figure  7b  assumes  a  daytime  steady-rate 
increase  of  0.1  per  hour  in  cd(t),  with  the  boundary  condi¬ 
tion  that  cd  =  cds  at  solar  noon,  where  cds  is  the  steady- 
state  value  of  cd,  as  shown  in  Figure  7a.  The  result  is  a 
better  match  of  servo  daytime  hQ  with  FLIP  hQ. 

The  assumption  of  a  daytime  steady  rate  for  cd  is  a 
type  of  servo  "fine  tuning."  The  USU  servo  model  makes 
this  tuning  optional.  Daytime  steady-rate  slope  values 
should  range  from  zero  (high  solar  radiation)  to  just  over 
0.1/hr  (low  solar  radiation).  FLIP  model  runs  indicate 
that  the  daytime  slope  depends  on  the  amount  of  solar 
radiation  available  for  ionization  production.  The  slope 
is  close  to  zero  for  low  geographic  latitudes  during  high 
solar  flux  conditions.  This  indicates  that  steady-state 
assumptions  are  valid  provided  there  is  sufficient  plasma 
production,  e.g.  in  the  summer  hemisphere. 

3.1.8  Neutral  Wind.  Once  the  vertically  applied 
plasma  drift  w  is  found  from  Equation  (18)  or  (19),  the  USU 
servo  model  calculates  the  meridional  neutral  wind  from 

-w 

U_  =  -  ( ?  8 ) 

p  cos  I  sin  I 
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where  Up  is  the  poleward  meridional  neutral  wind.  Equation 
(28)  assumes  the  total  plasma  drift  w  is  due  to  the  neutral 
wind.  In  general,  w  is  produced  by  both  the  neutral  wind 
and  electric  field.  If  the  electric  field  is  not  neglect¬ 
ed,  then  a  more  general  equation  for  the  upward  plasma 
drift  is 


w  =  Vj_pCos  I  -  Up  cos  I  sin  I  (29) 

where  Vj_p  is  the  ion  drift  velocity  component,  perpendicu¬ 
lar  to  the  geomagnetic  field,  in  the  poleward  direction, 
and  is  caused  by  an  eastward  electric  field  (Electric  field 
effects  are  explained  in  Section  2.1).  For  quiet  to  moder¬ 
ate  geomagnetic  conditions,  the  first  term  in  Equation  (29) 
can  generally  be  neglected,  leading  to  Equation  (28) . 
Equation  (28)  is  actually  preferred  for  many  ionospheric 
modeling  applications.  This  is  because  the  parameter  of 
interest  is  the  vertical  ion  drift,  regardless  of  how  that 
ion  drift  is  produced.  [Richards,  P.  G. ,  An  Improved  Algo¬ 
rithm  for  Determining  Neutral  Winds  from  the  Height  of  the 
F2-peak  Electron  Density,  Submitted  to  J.  Geophvs.  Res. . 
January,  1991]. 

3.2  Model  Program  Outline 

This  section  summarizes  the  equations,  parameters  and 
procedures  used  in  the  USU  servo  model  program.  The  fol¬ 
lowing  steps  are  presented  and  explained  in  sequence  and 
meant  as  a  guide  for  programming: 
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1.  Enter  the  location  parameters;  name,  geographic 
latitude,  longitude  and  dip  angle  (USU  servo  model  calcu¬ 
lates  an  optional  centered  dipole  dip  angle) . 

2.  Enter  the  date.  Solar  zenith  angles  (X)  are 
calculated  from  the  date,  and  geographic  coordinates,  which 
are  based  on  local  or  universal  time  [Blackadar,  1989]. 

3.  Enter  the  optional  value  for  the  daytime  steady- 
rate,  which  is  the  change  in  c^  per  hour,  or  else  a  zero 
slope  is  used  for  steady-state  approximation  (cd  remains 
constant) . 

4.  Input  MSIS  neutral  atmospheric  data  at  an  F2- 
region  reference  height  hr  (normally  300  km) .  This  may  be 
entered  as  an  input  data  file  or  MSIS  may  be  called  as  a 
subprogram. 

5.  Input  hmF2  data  file.  Here  the  IRI  may  be  used  as 
an  hmF2  data  source  from  a  file  or  subprogram. 

6.  Input  optional  ion  temperature  data  or  make  a 
ratio  assumption.  The  USU  servo  model  assumes  T^  =  l.lTn 
during  the  day,  =  1.05Tn  for  morning  or  evening  and  T^  = 
Tn  at  night. 

7.  Input  a  program  routine  to  find  when  the  day 
begins  and  ends  based  on  X.  This  step  is  only  needed  if  a 
daytime  steady-rate  was  entered  other  than  zero  to  change 
c^,  so  that  c^  is  not  constant  (step  3) . 

0.  The  program  determines  diurnal  periods  of  day, 
night,  morning  and  evening  based  on  X,  where  limits  are 


102°  and  82° . 
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9.  At  hr,  calculate  chemical  reaction  rates  k-^  and 
k2 ,  given  \n  Table  1. 

10.  At  hr,  calculate  loss  scale  height  factor  k  from 
Equation  ( £ ) . 

11.  At  hr,  and  geographic  latitude  calculate  gravity. 

12.  At  hr,  calculate  0+  scale  height  H. 

13.  At  hr,  calculate  linear  loss  coefficient,  intro¬ 
duced  in  Section  2.2  as  /3r  =  k^Oj]  +  k2[N2]. 

14.  At  hr,  calculate  ion-neutral  collision  frequency 
vr,  given  in  Table  2. 

15.  Use  Equation  (20)  to  calculate  the  boundary 
condition  balance  height  hQbc.  Use  cx  =  cd  =  1  and  cx  = 
acn  =  (2.82) (1.3).  For  twilight  periods,  linear  interpo¬ 
late  cx  between  cd  and  acn  with  time. 

16.  Use  Equation  (25)  to  calculate  hQ. 

17.  Use  Equation  (21)  to  calculate  the  servo  diurnal 
coefficients  cx,  which  are  now  "tuned"  with  an  hQ  input. 

18.  If  the  daytime  slope  from  step  3  is  not  equal  to 
zero,  then  recalculate  cd  =  cd(t).  From  step  17,  cx  =  cd  = 
cds,  where  cds  is  a  constant  steady-state  value.  Assume 
the  boundc.i  y  condition  that  cd(t)  =  cds  at  solar  noon, 
which  is  the  lowest  X  value.  Begin  to  recalculate  cd(t) 
when  the  day  period  first  begins  and  add  the  slope  value 
per  hour  so  as  to  pass  through  cds  at  solar  noon  and  end 
when  the  day  period  ends.  This  is  a  linear  change  in  the 
daytime  diurnal  coefficient  with  time. 

19.  If  the  daytime  slope  is  not  equal  to  zero,  then 
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use  Equation  (20)  to  recalculate  hQ  with  tuned  cx  values 
calculated  from  steps  17  and  18. 

20.  With  final  hQ  values,  calculate  the  ion-neutral 
collision  frequency  at  hQ,  where  vQ  =  vrexp(-zQr)  and 

zor  =  <ho  ~  hr>/H- 

21.  Use  Equation  (18)  or  (19)  to  calculate  the  verti¬ 
cal  plasma  drift  w. 

22.  Use  Equation  (28)  to  calculate  the  meridional 
neutral  wind  U. 

23.  At  this  point,  generate  an  output  file. 

Steps  1-23  are  generalized  and  presented  as  a  program 
guide  to  indicate  which  equations,  parameters  and  proce¬ 
dures  are  necessary  for  servo  calculations. 

3 . 3  Data  Preparation 

As  stated  in  Section  2.5,  the  USU  servo  model  requires 
MSIS  for  neutral  atmospheric  parameters  and  a  source  of 
hmF2  data.  MSIS  depends  on  location,  altitude,  day  of 
year,  solar  F10.7  cm  flux  and  geomagnetic  Ap,  all  given  in 
Appendices  A  and  B.  This  study  uses  three  sources  of  hmF2 
data  to  calculate  the  meridional  neutral  wind.  The  first 
source  is  from  the  IRI,  which  generates  monthly  averages  of 
hmF2  data  for  magnetically  quiet  conditions  in  the  non- 
auroral  ionosphere.  The  second  source  is  provided  by  20 
ionosonde  stations,  which  were  operating  during  a  world  day 
campaign  from  12-16  January  1988  (GISMOS-88).  The  third 


source  of  hmF2  data  comes  from  incoherent  scatter  radar 
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measurements,  from  the  Millstone  Hill  Haystack  Observatory, 
Massachusetts,  during  the  ETS  period  of  study,  from  17-24 
September  1984. 

3.3.1  IRI  hmF2  Data.  The  IRI  is  used  to  calculate  two 
sets  of  data  to  show  two  types  of  application.  Appendices 
A  and  B  have  the  input  parameters  for  location,  12-month- 
mean  sunspot  number  and  month.  The  first  example  shows  a 
neutral  wind  comparison  between  low  and  high  solar  activity 
during  geomagnetically  quiet  periods,  for  two  magnetically 
conjugate  locations.  The  two  locations  used  in  the  solar 
cycle  application  are  Petropavlovsk,  USSR  and  Canberra, 
Australia.  The  second  example  from  the  IRI  is  used  to 
calculate  neutral  winds  along  a  meridional  chain  of  iono- 
sonde  stations  in  the  northern  and  southern  hemispheres. 
This  second  IRI  data  set  was  picked  to  correlate  with  data 
from  20  ionosonde  stations  that  made  measurements  during 
the  GISMOS-88  world-day  campaign. 

3.3.2  Ionosonde  hmF2  Data.  F2-peak  heights  are  deter¬ 
mined  using  ionosonde  data  from  an  empirical  relationship 
developed  by  Bradley  and  Dudeney  [1973],  as  discussed  in 
Section  2.5.  The  three  necessary  elements  of  data,  namely 
foE,  f oF2  and  M(3000)F2  are  from  the  scientific  archives  in 
the  World  Data  Center  at  Boulder,  Colorado.  As  stated, 
there  are  problems  in  measuring  foE  at  night,  so  that  the 
E-layer  critical  frequency  is  calculated  from  the  IRI  model 
E  layer  density  NmE,  where  foE  =  9  x  lO'^NmE)1.  Values  of 
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foE  are  used  in  a  ratio  of  foF2  to  foE,  to  compare  the 
relative  importance  of  the  F  and  E-layer  ionization  in  the 
calculation  of  hmF2 .  The  variance  between  the  IRI  and 
measured  values  of  foE  is  normally  less  than  10%  and  does 
not  significantly  impact  the  neutral  wind  calculation. 
Appendix  A  lists  the  20  ionosonde  locations  and  Appendix  B 
shows  the  GISMOS-88  ionospheric  parameters. 

3.3.3  Incoherent  Scatter  Radar  Data.  For  this  study, 
Millstone  Hill  incoherent  scatter  radar  data  was  graciously 
provided  by  Dr.  Michael  J.  Buonsanto  of  the  Millstone  Hill 
Haystack  Observatory.  The  radar  data  covers  the  ETS  period 
of  September  17-24,  1984.  The  radar  was  operated  in  a  mode 
that  used  steerable  north-south  elevation  scans  and  a  fixed 
elevation  measurement.  The  method  of  radar  data  collection 
and  neutral  wind  derivation  is  fully  discussed  by  Buonsanto 
et  al.  [1989].  Included  in  the  radar  data  is  hmF2 ,  Ti,  Te, 
radar-derived  neutral  winds  and  statistical  uncertainties 
in  the  radar  winds.  Buonsanto  et  al.  [1989]  used  this  same 
ETS  period  to  make  neutral  wind  comparisons  of  radar-de¬ 
rived  winds  with  Buonsanto  et  al.  [1989]  servo  model  winds 
and  Miller  et  al.  [1986]  FLIP  model  winds.  The  ionospheric 
input  parameters  for  the  ETS  period  are  presented  in  Appen¬ 


dix  B. 
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CHAPTER  IV 

RESULTS  AND  DISCUSSION 

4.1  Winds  from  IRI  hmF2 

Changes  in  mid-latitude  neutral  wind  speeds  depend  on 
changes  in  momentum  forcing  terms  in  Equation  (1).  The 
pressure  gradient  force  is  the  main  driver  of  the  neutral 
winds  and  is  caused  by  temperature  gradients  created  by  UV 
and  EUV  radiation  heating.  An  increase  in  solar  activity 
directly  increases  radiation  levels  and  intensity,  thereby 
increasing  horizontal  pressure  gradients.  Increases  in 
solar  activity  also  leads  to  an  increase  in  charged  parti¬ 
cle  population,  which  leads  to  an  increase  in  ion  drag. 

The  increase  in  pressure  gradient  force  and  increase  in  ion 
drag  are  counterproductive  in  the  movement  of  neutral  air. 
In  general,  as  solar  activity  increases,  neutral  winds  are 
expected  to  decrease  because  ion  drag  increases  more  rapid¬ 
ly  than  the  pressure  gradient  force  [Babcock  and  Evans, 
1979;  Rishbeth,  1972]. 

Figure  8  shows  meridional  neutral  winds  calculated  for 
two  magnetically  conjugate  locations  during  January,  1981 
and  1985.  Solar  activity  was  high  in  1981  and  low  in  1985. 
Values  of  hmF2  were  modeled  using  the  IRI.  The  dashed 
lines  represent  neutral  winds  for  high  solar  activity, 
while  the  solid  lines  depict  low  solar  activity.  Iono¬ 
spheric  parameters  for  both  periods  are  given  in  Appendix 
B.  Figure  8  shows  that  nighttime  neutral  wind  speeds  are 
significantly  smaller  during  high  solar  activity.  When  a 
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Fig.  8.  Meridional  neutral  winds  for  two  magnetically 
conjugate  locations  during  low  (January,  1985)  and  high 
solar  activity  (January,  1981)  shown  as  solid  and  dashed 
lines  respectively.  Winds  are  calculated  by  the  USU  servo 
model  using  IRI  F2-peak  heights. 
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neutral  particle  strikes  an  ion  moving  at  a  different 
velocity,  some  of  the  neutral's  momentum  is  transferred  to 
the  ion  and  the  ion  velocity  is  changed.  The  momentum  lost 
to  the  neutral  particle  is  an  ion-drag  force  on  the  neu¬ 
tral.  Ion  drag  increases  with  increasing  ionized  particle 
population.  Figure  8  shows  that  increased  solar  activity 
decreases  nighttime  equatorward  winds  by  as  much  as  50  m  s' 
at  Canberra  and  more  than  75  m  s'1  at  Petropavlovsk.  This 
decrease  is  attributed  to  a  larger  nighttime  ion  drag  at 
high  solar  activity.  Figure  8  also  shows  a  marked  increase 
in  the  poleward  duration  of  the  neutral  winds  during  high 
solar  activity,  but  little  decrease  in  wind  speed  from  low 
to  high  solar  activity. 

Figure  8  shows  little  change  in  daytime  wind  speeds 
when  comparing  wind  patterns  at  low  and  high  solar  activi¬ 
ty.  Although  ion  production  increases  daytime  ion  drag, 
pressure  gradients  also  increase  during  the  day  and  at  high 
solar  activity.  In  this  case,  Figure  8  indicates  that  an 
increase  in  daytime  pressure  gradients  at  high  solar  activ¬ 
ity  tends  to  counterbalance  the  increased  ion  drag  and 
therefore  little  change  is  observed  in  the  daytime  wind 
speeds. 

Figure  9  shows  another  example  of  IRI  produced  winds. 
In  this  example,  a  meridional  chain  of  10  ionosonde  sta¬ 
tions  show  latitudinal  changes  in  neutral  wind  patterns. 
These  wind  patterns  are  calculated  for  the  GISMOS-88  period 
to  allow  future  comparison  with  ionosonde  produced  winds. 
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Fig.  9.  Neutral  winds  along  a  meridional  string  of  iono- 
sonde  stations  calculated  by  the  USU  servo  model  from  IRI 
hmF2  data  for  January,  1988.  Positive  winds  are  poleward 
Open  and  shaded  circles  mark  local  noon  and  midnight. 
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As  before,  the  diurnal  patterns  of  daytime  poleward  winds 
and  nighttime  equatorward  winds  are  quite  plain.  Of  note 
is  the  evidence  of  a  much  stronger  diurnal  tide  and  greater 
magnitude  of  meridional  neutral  wind  speeds  at  the  higher 
latitudes.  This  is  probably  due  to  the  much  larger  diurnal 
variation  of  atmospheric  heating  at  higher  latitudes. 

Miller  et  al.  [1989]  made  a  similar  comparison  of  IRI  and 
ionosonde  winds.  Comparative  results  from  the  Miller  et 
al.  [1989]  study  shows  that  IRI  winds  are  nearly  the  same 
as  ionosonde  winds  derived  from  median  ionosonde  hmF2 
values  over  a  15-day  period. 

4.2  Winds  from  Ionosonde  ..  ^2 

The  World  Day  campaign  from  12-16  January,  1988,  also 
called  the  GISMOS-88  campaign,  is  used  for  a  basis  of 
study.  For  this  application,  20  ionosonde  stations  along 
an  approximate  meridional  longitude  of  135°  E  geographic 
are  used  to  show  a  latitudinal  variation  in  the  neutral 
winds.  A  moderately  strong  geomagnetic  storm  occurred 
about  half  way  through  the  GISMOS-88  campaign.  For  com¬ 
parison,  the  neutral  wind  plots  are  broken  into  prestorm, 
storm  and  post-storm  periods  of  12-13,  14-15  and  16  January 
respectively.  Table  3  shows  the  change  in  daily  Ap  from  9- 
17  January.  The  changes  in  Ap  are  used  to  explain  changes 
in  neutral  wind  circulation  patterns. 

Figure  10  shows  the  meridional  neutral  wind  patterns 
for  two  days  prior  to  the  storm  event.  Dashed  lines  are 
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Fig.  10.  Prestorm  neutral  winds  (positive  poleward)  for  c. 
meridional  string  of  20  ionosonde  stations  on  January  12 
and  13,  19P8  (GISMOS-88 ) .  Winds  are  calculated  by  the  USU 
servo  model. 
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Fig.  10.  (continued) 
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Fig.  10.  (continued) 
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Fig.  10.  (continued) 


TABLE  3.  Geomagnetic  Ap  from  9-37  January,  1988 

Date : 

Ap: 

used  to  join  the  wind  pattern  at  points  where  hourly  obser¬ 
vations  are  missing.  Although  the  12th  and  13th  are  rela¬ 
tively  quiet  compared  to  the  storm  days,  there  are  still 
distinguishable  differences  between  the  two  days.  The 
differences  are  probably  attributable  to  differences  in 
geomagnetic  conditions.  On  the  12th,  the  wind  pattern  is 
more  erratic  and  nighttime  equatorward  winds  are  greater  at 
most  locations  in  the  northern  and  southern  hemispheres. 

The  stronger  equatorward  winds  are  most  prevalent  above 
3 5  ° -4 0 ° ,  suggesting  that  increases  in  auroral  heating 
processes  may  have  enhanced  equatorward  pressure  gradients. 
There  may  also  have  been  enhanced  equatorward  winds  at  the 
higher  latitudes  due  to  increased  ion-neutral  momentum 
transfer  from  the  polar  cap  convective  cells. 

In  Figure  11,  much  of  the  hourly  observations  are 
unfortunately  not  available.  This  is  probably  due  to 
ionospheric  short  wave  fadeout  from  the  storm  event  [Har¬ 
greaves,  1979;  Ratcliffe,  1972].  Most  of  the  fadeout  is 
probably  caused  by  enhanced  D-layer  ionization.  This 
condition  normally  affects  higher  latitudes,  but  missing 
data  is  evident  throughout  almost  the  entire  ionosonde 
network.  Again  dashed  lines  are  used  to  close  the  gap 
between  known  observations. 


10  11  12  13  14  15  16  17 

4  11  21  7  48  63  5  7 


YAKUTSK  62.  0  N,  129.  6  E 


PETROPAVLOVSK  53.0  N,  15  8.6  E 


0  6  12  18  24  30  36  42 


UT  (hours) 


Fig.  11.  Neutral  winds  during  GISMOS-88  storm  on  January 
14  and  15,  1988  for  a  meridional  string  of  20  ionosonde 
stations.  Winds  calculated  by  USU  servo  model. 
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Fig.  11.  (continued) 
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The  most  striking  feature  in  Figure  11  is  a  change  in 
the  daytime  neutral  wind  pattern  on  the  15th.  At  most 
locations,  daytime  winds  are  either  close  to  zero  or  equa- 
torward.  Much  of  the  data  is  missing  on  the  night  of  the 
storm,  between  the  14th  and  15th.  However,  there  is  still 
indications  of  minor  equatorward  surges,  especially  at 
higher  latitudes.  The  magnitudes  of  the  equatorward  w;nd 
surges  are  not  as  strong  as  expected  when  compared  with  the 
two  relatively  quiet  days  of  the  12th  and  13th.  This  may 
be  explained  by  considering  dynamo  electric  fields. 

A  strong  equatorward  wind  surge  during  the  storm  would 
cause  a  westward  electric  field  during  the  build-up  of  the 
surge.  This  electric  field  drives  the  plasma  downward, 
perpendicular  to  the  field  lines.  Using  hmF2  to  calculate 
meridional  neutral  winds  during  a  storm  will  therefore  give 
a  wind  speed  less  than  it  really  is.  Blanc  and  Richmond 
[1980]  and  Foster  and  Aarons  [1988]  described  this  type  of 
disturbance  dynamo  generated  by  equatorward  neutral  winds 
driven  by  enhanced  auroral  heating  at  high  latitudes.  In 
Section  4.3,  comparisons  of  neutral  wind  techniques  at 
Millstone  Hill  show  that  strong  westward  electric  fields  at 
the  time  of  storm  occurrence  are  responsible  for  most  of 
the  neutral  wind  correction  when  using  the  hmF2  method  in 
calculating  neutral  winds. 

Figure  12  shows  the  post-storm  wind  patterns.  For  the 
most  part,  neutral  winds  have  returned  to  normal  diurnal 
patterns.  Overall  poleward  winds  on  the  16th  are  a  little 
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more  enhanced  than  the  prestorm  period.  An  increase  in  the 
uncertainties  in  MSIS-86  neutral  densities  and  temperatures 
after  the  storm  event  may  account  for  the  differences  in 
pre  and  post-storm  wind  patterns,  rather  than  a  post-storm 
dynamic  process. 

The  IRI  neutral  winds  in  Figure  9  presents  a  climato¬ 
logical  background  for  comparison  with  Figures  10,  11  and 
12.  The  general  patterns  in  the  IRI  winds  are  reflected  on 
the  patterns  of  the  ionosonde  winds.  Fluctuations  depart¬ 
ing  from  the  climatological  patterns  may  be  caused  by 
" ionospheric  weather . " 

Figure  13  shows  the  ionosonde  winds  for  the  prestorm 
period  (dashed  line)  at  Petropavlovsk  superimposed  on  the 
IRI  wind  pattern  (solid  line) .  Day  13  is  the  more  quiet  of 
the  two  days  and  the  ionosonde  pattern  tends  to  follow  the 
IRI  pattern  somewhat  closer.  Changes  in  neutral-wind  or 
electric-field  induced  plasma  motions  may  be  responsible 
for  the  differences  in  the  two  patterns. 

Figure  13  shows  how  IRI  wind  patterns  approximate 
daily  observed  ionosonde  wind  patterns.  Figure  13  also 
shows  that  IRI  wind  patterns  may  be  used  as  a  climatologi¬ 
cal  basis  in  the  study  of  transcient  neutral-wind  and 
electric-field  features. 

4.3  Comparing  Winds  with  Radar  Data 

The  equinox  transition  study  during  September  17-24, 
1984,  provides  a  period  of  study  in  which  to  compare  the 
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Fig.  13.  Weather  vs.  Climatology:  Neutral  winds  from  the 
USU  servo  model  for  January  12  and  13,  1988  at  Petropav- 
lovsk,  during  GISMOS-88  prestorm  period.  Solid  lines  show 
IRI  winds  and  dashed  lines  show  ionosonde  winds. 
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USU  servo  model  with  other  wind  calculation  techniques, 
during  a  variety  of  geomagnetic  conditions.  The  first  part 
of  the  ETS  period  was  geomagnetically  quiet.  A  small  storm 
(Ap  =  36)  occurred  on  19  September,  followed  by  a  three  day 
recovery  period  and  then  a  major  storm  (Ap  =  112)  on  23 
September . 

Buonsanto  et  al.  [1989]  used  this  same  period  of  study 
and  made  neutral  wind  comparisons  between  radar,  the  Miller 
et  al.  [1986]  FLIP  method  and  servo  method  adapted  by 
Bounsanto  et  al.  [1989].  To  distinguish  between  the  wind 
techniques,  we  will  refer  to  the  methods  of  calculation  as 
FLIP  winds,  USU  servo  winds,  servo  winds  and  radar  winds. 

The  FLIP  and  servo  winds  are  similar  during  this  ETS 
period  of  study,  but  the  servo  winds  are  consistently  more 
negative  (southward) .  The  servo  model  is  tuned  for  better 
agreement  with  the  radar  winds,  which  accounts  for  less 
than  half  this  difference  between  FLIP  and  servo  winds, 
since  the  FLIP  winds  are  not  radar  tuned.  Differences 
between  the  servo  and  FLIP  winds  are  greatest  during  geo¬ 
magnetically  active  times.  This  difference  in  winds  be¬ 
tween  the  two  techniques  partially  prompted  this  investiga¬ 
tion. 

Figure  14  shov/s  neutral  winds  above  Millstone  Hill 
during  18,  19  and  20  September  calculated  by  the  FLIP,  USU 
servo  and  servo  models.  Local  midnight  at  Millstone  Hill 
corresponds  to  0446  UT.  Daily  Ap  history  is  input  into 
MSIS-86  and  1.7  is  used  to  increase  the  collision  cross 
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Fig.  14.  Neutral  wind  comparison  above  Millstone  Kill 
during  ETS  period.  All  calculations  include  d?>  Lly  Ap 
history  and  1.7  increase  to  0+->0  collisions.  Solid  line  is 
FLIP  winds,  short-dashed  line  is  USU  servo  winds  and  long- 
dashed  line  is  tuned  servo  model  winds. 
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section  of  0+-"0  for  neutral  wind  calculations.  The  solid 
line  depicts  FLIP  winds,  short  dashes  are  for  USU  servo 
winds  and  long  dashes  are  for  servo  winds.  Measured  is 
used  in  USU  servo,  while  measured  and  Te  are  used  to 
calculate  servo  winds.  The  FLIP  model  internally  calcu¬ 
lates  its  own  plasma  temperatures.  The  USU  servo  and  servo 
models  assume  daytime  steady-state  conditions,  where  cd 
remains  constant,  as  discussed  in  Section  3.1.7.  The  servo 
model  is  tuned  with  cd  =  1  and  cn  =  1.3  as  in  Buonsanto  et 
al.  [1989].  The  USU  servo  model  requires  no  tuning  because 
it  uses  a  balance  height  adjustment  technique  described  in 
Section  3.1.5.  We  can  input  the  USU  servo  adjusted  hQ 
values  into  Equation  (21)  or  (22)  and  calculate  equivalent 
tuning  coefficients  to  compare  with  the  Buonsanto  et  al. 
[1989]  values.  These  USU  servo  equivalent  tuning  coeffi¬ 
cients  are  cd  =  1.05  and  cn  =  1.41. 

Results  from  Figure  14  show  that  the  USU  servo  winds 
are  in  closer  agreement  to  the  FLIP  winds  than  are  the 
servo  winds.  The  USU  servo  winds  are  5  to  20  m  s'1  more 
southward  than  the  FLIP  model  winds,  while  the  servo  winds 
are  between  20  to  50  m  s'1  more  southward  on  the  quiet  day. 
On  the  active  day  of  the  19th  and  continuing  into  the  post¬ 
storm  day  of  the  20th,  the  separation  between  FLIP  and 
servo  winds  grow  to  about  50  to  70  m  s'1.  The  USU  servo 
winds  maintain  about  the  same  separation  throughout  the 
entire  period. 

There  are  two  main  reasons  why  the  USU  servo  winds  are 
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closer  to  the  FLIP  winds.  The  lesser  reason  is  that  the 
servo  winds  are  empirically  tuned  to  create  greater  south¬ 
ward  wind  speeds,  which  are  in  better  agreement  with  the 
radar  winds.  The  second  and  more  physical  reason,  stems 
from  the  replacement  of  the  peak  vertical  plasma  velocity 
Vzm  in  Equation  (11).  In  the  USU  servo  derivation,  Equa¬ 
tion  (7)  is  V2m  =  w  -  (g/vm)sin2I.  In  the  RGW  servo  model, 
which  is  used  by  Buonsanto  et  al.  [1989],  V2m  =  w  - 
(Dpm/2H)sin2 I. 

Two  improvements  come  about  by  the  use  of  Equation 
(7) .  One  is  that  vm  depends  on  T^ ,  while  Dpm  depends  on 
both  T^  and  Te,  but  the  second  and  more  important  improve¬ 
ment  is  that  Equation  (7)  does  not  assume  diffusive  equi¬ 
librium.  At  diffusive  equilibrium,  Tn  =  T^  =  Te  and  Hp  may 
be  replaced  by  2H,  as  discussed  in  Section  2.4.  The  re¬ 
placement  of  Vzm  in  Equation  (11)  affects  both  plasma  drift 
w  from  Equation  (18)  or  (19)  and  the  balance  height  hQ  from 
Equation  (20) . 

To  explain  the  large  differences  between  the  FLIP  and 
servo  winds  during  geomagnetically  active  periods,  consider 
Equation  (19),  rewritten  here  and  followed  by  an  RGW  servo 
model  equivalent 


g  sin2 1 

w  =  -  zm(k  +  1) 


,  Dpp  sin2 I 


(30) 


W 


2H 


zm(k  +  1) 


(31) 
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At  diffusive  equilibrium,  both  Equations  (30)  and  (31)  are 
equal.  During  periods  of  increased  Te  and  Tj_,  diffusive 
equilibrium  no  longer  applies  and  the  two  equations  di¬ 
verge.  As  plasma  temperatures  increase,  w  decreases  in 
Equation  (30)  and  increases  in  Equation  (31) .  Physically, 
Equation  (30)  is  more  correct  because  collision  frequency 
increases  with  increasing  temperature.  As  collisions 
between  the  ions  and  neutrals  increase,  the  neutral  wind 
becomes  more  efficient  in  moving  the  plasma  along  the  field 
lines  and  w  decreases.  Since  the  scale  height  H  in  the 
denominator  of  Equation  (31)  is  of  atomic  oxygen,  then  the 
factor  2  may  be  replaced  by  (T^  +  Tg)/Tn  to  become  the 
plasma  scale  height  Hp.  An  increase  in  plasma  temperature 
wou 1 d  then  decrease  w  in  Equation  (31)  and  equal  Equation 
(30)  . 

Figure  15  shows  the  neutral  and  ion  temperatures  for 
18,  19  and  20  September  and  Figure  16  shows  Te  for  the  same 
period.  and  Te  are  measured  with  radar  at  300  km,  while 

Tn  is  from  MSIS-86  at  300  km.  Plasma  temperature  increases 
during  and  after  the  storm  event  on  the  19th  and  servo 
winds  in  Figure  14  have  grown  consistently  more  southward. 

How  much  does  an  increase  in  affect  the  USU  servo 
model  winds?  Figure  17  shows  USU  servo  model  winds  calcu¬ 
lated  from  assumed  ratios  of  Tj/Tn  =1,  1.5  and  2.5  on  18 
September,  depicted  by  solid,  short-dashed  and  long-dashed 
lines  respectively.  The  hmF2  radar  data  used  for  the 
calculations  shown  in  Figure  17  are  from  repetitive  eleva- 
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Fig.  15.  Neutral  and  ion  temperatures  at  300  km  above 
Millstone  Hill,  during  ETS  period.  T„  is  solid  line  from 
MSIS-86  and  T j  is  dashed  line  measured  by  IS  radar. 
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Fig.  17.  USU  servo  neutral  wind  sensitivity  to  changes  in 
ion  temperature  during  Millstone  Hill  ETS  period.  Winds 
are  calculated  for  different  temperature  ratios.  Solid 
line  is  for  Tj/Tn  =  1,  short-dashed  line  is  for 
Tj/Tn  =  1.5,  and  long-dashed  line  is  for  Tj/Tn  =  2.5. 
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tion  scans  and  have  a  higher  resolution  than  the  data  used 
in  Figures  14  or  18.  Figure  17  shows  that  higher  values  of 
T^  tend  to  decrease  the  neutral  wind,  due  to  an  increase  in 
collision  frequency. 

Considering  the  large  differences  in  used  in  the 
neutral  wind  calculations  of  Figure  17,  ion  temperature  is 
not  a  major  contributor  in  the  uncertainty  of  the  USU  servo 
model.  We  can  therefore  assume  the  uncertainty  in  the  USU 
servo  model  from  uncertainties  in  is  small  compared  with 
uncertainties  from  electric  field  effects  during  active 
geomagnetic  periods.  This  is  fortunate,  since  assumed  ion 
temperatures  are  used  in  calculating  neutral  winds  from 
ionosonde  hmF2  data. 

The  effect  of  electric  fields  on  plasma  drift  is 
significant  during  a  geomagnetic  storm.  The  perpendicular 
ion  drift  velocity  term  in  Equation  (29)  must  be  included 
in  calculations  of  neutral  winds  obtained  by  an  hmF2  method 
[Buonsanto  et  al.,  1989].  Figure  18  compares  winds  calcu¬ 
lated  by  the  USU  servo  model  and  incoherent  scatter  radar 
at  Millstone  Hill  from  18-23  September.  Radar  winds  are 
calculated  at  300  km  by  the  Salah  and  Holt  [1974]  method, 
with  some  modifications  discussed  by  Buonsanto  et  al. 
[1989].  The  1.7  increase  in  0+-»0  collision  cross  section 
is  used  in  both  the  USU  servo  and  radar  winds.  The  three 
hourly  ap  history  is  used  in  the  radar  calculation  and 
daily  Ap  history  is  used  for  the  USU  servo  calculation. 

The  use  of  ap  history  generally  gives  stronger  equatorward 
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winds,  especially  at  night  and  during  the  geomagnetic  storm 
periods  [Buonsanto  et  al.,  1989].  Therefore,  a  slight 
equatorward  bias  in  radar  derived  winds  are  evident  when 
compared  with  the  USU  servo  winds.  Error  bars  indicate  the 
statistical  uncertainties  in  the  radar  wind  and  are  used  to 
gauge  the  USU  servo  wind  approximation.  Uncertainties  in 
the  USU  servo  model  are  not  considered  in  this  comparison, 
but  as  a  minimum  should  approximate  the  same  error  as 
reflected  by  the  radar. 

Except  for  the  storm  times,  USU  servo  winds  generally 
fall  within  the  radar  wind  error  bars.  This  approximate 
agreement  indicates  that  the  balance  height  adjustment 
technique  works  and  tuning  the  USU  servo  model  is  unneces¬ 
sary. 

Figure  18  shows  large  differences  in  the  USU  servo  and 
radar  winds  during  both  storms.  The  mild  storm  on  the  19th 
lasted  for  about  6  hours,  while  the  strong  storm  on  the 
23rd  lasted  for  about  15  hours.  Ion  drift  velocity  meas¬ 
urements  calculated  by  Buonsanto  et  al.  [1989]  gives  elec¬ 
tric  field  corrections  of  nearly  200  m  s'1  for  the  storms. 
This  electric  field  correction,  combined  with  increased 
uncertainties  of  about  100  m  s*1  for  both  radar  and  USU 
servo  model  calculations  during  geomagnetic  storm  events 
[Buonsanto,  et  al.,  1989,  1990],  explains  the  observed 
large  difference  in  the  USU  servo  and  radar  winds. 
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CHAPTER  V 

CONCLUSIONS  AND  RECOMMENDATIONS 

At  mid-latitude  thermospheric  heights,  solar  UV  and 
EUV  radiation  heating  contributes  most  of  the  energy  to 
drive  the  neutral  winds.  This  level  of  solar  radiation 
increases  directly  with  solar  activity.  As  solar  activity 
increases,  an  increase  in  temperature  variation  will  in¬ 
crease  the  horizontal  pressure  gradient,  which  is  the  main 
driver  of  the  mid-latitude  thermospheric  winds. 

An  increase  in  solar  radiation  not  only  produces  heat, 
but  also  increases  the  charged  particle  population.  Colli¬ 
sions  between  ions  and  neutral  particles  become  more  fre¬ 
quent,  so  ion  drag  increases  with  increasing  solar  activi¬ 
ty.  The  pressure  gradient  and  ion  drag  are  counterproduc¬ 
tive  in  neutral-air  motion  and  lead  to  a  change  in  the 
balance  of  forces.  The  study  of  neutral  wind  force  balance 
relationships  from  solar  cycle  variations  leads  to  a  better 
understanding  of  mid-latitude  thermospheric  circulation 
patterns.  The  use  of  IRI  and  ionosonde  hmF2  data  in  the 
USU  servo  model  provides  a  valuable  means  of  calculating 
global  scale  variations  in  neutral  wind  patterns. 

The  USU  servo  model  uses  a  balance  height  adjustment 
technique  to  compensate  for  changes  in  chemical  and  diffu¬ 
sive  loss  of  ions  at  different  mid-latitude  locations  and 
ionospheric  conditions.  This  balance  height  adjustment 
technique  physically  replaces  servo  model  empirical  tuning 
by  radar  [Buonsanto  et  al.,  1989;  1990]  and  greatly  in- 
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creases  the  flexibility  in  the  servo  calculation  method. 

The  USU  servo  model  is  physically  more  realistic  than 
previous  servo  models.  In  the  USU  servo  model,  the  verti¬ 
cal  ion  velocity  V2m  is  calculated  by  w  -  (g/vm)sinJl 
instead  of  w  -  (Dpm/2H) sin2 I ,  as  in  the  theoretical  RGW 
servo  model.  The  2H  in  the  RGW  servo  model  assumes  diffu¬ 
sive  equilibrium  and  is  not  physically  realistic.  Electron 
and  ion  temperatures  are  used  to  calculate  Dpm,  but  only 
ion  temperature  is  required  to  calculate  vm,  making  this 
Vzm  replacement  a  more  practical  one. 

The  IRI  empirically  models  monthly  averages  of  hmF2 
for  magnetically  quiet  conditions  in  the  non-auroral  iono¬ 
sphere.  Since  the  IRI  gives  average  values  of  hmF2  for 
quiet  conditions,  its  use  is  ideal  in  che  calculation  of 
neutral  winds  from  hmF2  methods.  Random  errors  and  elec¬ 
tric  field  variability,  as  compared  to  solar  driven  diurnal 
forces,  are  averaged  out  of  the  IRI  hmF2  values.  This 
allows  calculation  of  global  climatological,  mid-latitude 
neutral  winds  for  different  solar  conditions. 

The  USU  servo  model  provides  a  fast  and  comparatively 
accurate  means  of  calculating  hmF2  neutral  winds  for  any 
mid-latitude  location,  without  regard  to  tuning.  I  recom¬ 
mend  that  a  servo  subprogram  using  the  equations  and  proce¬ 
dures  as  described  in  this  study,  be  implemented  as  part  of 
the  IRI  program. 

The  generation  of  neutral  winds  from  hmF2  proves  to  be 
a  valuable  technique.  The  major  uncertainty  is  the  result 
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of  induced  plasma  motions  from  electric  fields.  Under 
quiet  to  moderate  geomagnetic  conditions,  uncertainties 
from  electric  fields  are  usually  less  than  the  statistical 
uncertainties  in  wind  calculations.  For  active  geomagnetic 
periods,  one  should  consider  larger  uncertainties  produced 
by  increased  electric  fields. 

Figure  18  gives  an  indication  of  uncertainty  magni¬ 
tudes  and  duration  in  USU  servo  and  radar-derived  winds 
during  mild  and  strong  storm  events  at  Millstone  Hill.  The 
duration  of  large  uncertainties  from  electric  field  effects 
for  the  two  ETS  storms  is  about  6  hours  for  the  mild  storm 
and  about  15  hours  for  the  strong  storm.  Buonsanto  et  al. 
[1989]  calculated  electric  field  wind  corrections,  which 
compensate  for  most  of  the  uncertainty  between  the  radar 
and  USU  servo  model  winds  in  Figure  18.  If  electric  field 
effects  are  not  included  in  hmF2  neutral  wind  calculations 
during  geomagnetically  active  periods,  then  neutral  winds 
will  include  errors  exceeding  model  uncertainties. 

A  recent  paper  by  Sica  and  Schunk  [1990]  showed  that 
the  ionosphere/plasmasphere  vertical  exchange  flux  is 
important  in  modeling  nighttime  hmF2  values  if  only  the 
peak  ion  density  is  known.  The  topside  shape  factor  a,  as 
defined  by  Equation  (10)  and  discussed  in  Section  3.1.6, 
also  affects  F2-layer  height  and  may  be  related  to  The 

topside  electron  density  profile  can  be  measured  by  inco¬ 
herent  scatter  radar  and  used  to  compute  a  and  perhaps 
approximate  Changes  in  will  affect  the  nighttime 
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decay  rate  and  a  should  change  as  a  result.  To  determine  a 
relationship  between  the  two  parameters,  Equation  (21)  or 
(22)  may  be  used  to  calculate  nighttime  cx  values  from  the 
FLIP  model  for  varying  values.  The  nighttime  value  of 
cx  is  more  sensitive  to  changes  in  decay  rates  than  hmF2 , 
making  cx  a  more  desirable  quantity  in  the  study  of  verti¬ 
cal  plasma  flux  effects. 

The  FLIP  model  indicates  that  the  daytime  ion  produc¬ 
tion  should  be  approximated  by  a  steady-rate,  rather  than  a 
steady-state  assumption,  especially  during  periods  of  low 
solar  ionization  (see  Section  3.1.7).  This  allows  for  a 
daytime  continuity  equation  of  the  form  qm  =  cd(t)/3mN_, 
where  cd(t)  changes  at  a  steady  daytime  rate.  How  fast 
this  change  occurs  is  a  function  of  solar  ion  production, 
which  depends  on  geographic  latitude  (or  solar  zenith 
angle)  and  solar  radiation  flux.  The  FLIP  model  indicates 
that  cd(t)  approaches  a  constant  value  as  daytime  ion 
production  increases,  implying  steady-state  conditions. 

The  USU  servo  model  gives  the  option  of  inputting  a  steady- 
rate  in  cd(t)  per  hour.  I  recommend  the  development  of  a 
cd(t)  linear  change  function  based  on  solar  zenith  angle 
and  solar  flux  intensity  (perhaps  a  Chapman  type  function) . 
The  daytime  slope  function  may  then  be  included  in  the  USU 
servo  model  for  a  physically  better  daytime  output. 
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Appendix  A.  General  Information  on  Station  Locations 


TABLE  4.  Ionosonde/Radar  Station  Parameters 


Name 

Geoaraohic  Location 

Dip  Anal' 

YAKUTSK 

62.00° 

N 

129.60° 

E 

76.1° 

MAGADAN 

60.00° 

N 

151.00° 

E 

71.1° 

PETRO PAVLOVS K 

53.00° 

N 

158.60° 

E 

64.3° 

MANCHOULI 

49.60° 

N 

117.40° 

E 

67.8° 

KHABAROVSK 

48.50° 

N 

135.10° 

E 

63.6° 

WAKKANAI 

45.40° 

N 

141.70° 

E 

59.3° 

AKITA 

39.70° 

N 

140.10° 

E 

53.4° 

BEIJING 

39.90° 

N 

116.40° 

E 

57.7° 

TOKYO 

35.70° 

N 

139.50° 

E 

48.7° 

YAMAGAWA 

31.20° 

N 

130.60° 

E 

44.1° 

OKINAWA 

26.30° 

N 

127.80° 

E 

36.8° 

CHUNG-LI 

24.91° 

N 

121.24° 

E 

35.1° 

CANTON 

23.10° 

N 

113.30° 

E 

32.3° 

VANIMO 

2.70° 

S 

141.30° 

E 

-22.5° 

DARWIN 

12.45° 

S 

130.95° 

E 

-40.7° 

TOWNSVILLE 

19.63° 

S 

146.85° 

E 

-49.3° 

MUNDARING 

31.98° 

s 

116.22° 

E 

-67.1° 

CANBERRA 

35.32° 

s 

149.00° 

E 

-66.4° 

SALISBURY  (AUST.) 

34.70° 

s 

138.60° 

E 

-67.2° 

HOBART 

42.92° 

s 

147.32° 

E 

-72.9° 

Millstone  Hill 

42.50° 

N 

288.50° 

E 

71.5° 

Magnetic  Dip  Angles  from  International  Geomagnetic 
Reference  Field  (IGRF) 
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Geomagnetic  and  Solar  Data 


TABLE  5.  Geomagnetic  and  Solar  Data 


1985 

1981 


1984 

1984 

1984 

1984 

1984 

1984 

1984 

1984 


F10. 7 


100.9 

101.7 

107.5 
108.1 

113.7 
112.4 

121.8 


72.6 

161.6 


73.4 

74.6 

73.8 

74.6 

74.1 

75.1 

75.9 

76.1 


104.6 

104.6 

104.6 

104.6 

104.6 

104.6 

104.6 


73 . 1 
205.1 


78.9 

78.9 

78.9 

78.9 

78.9 

78.9 

78.9 

78.9 


A  Ad 

— lz 

4 

32.8 

11 

32.8 

21 

32 . 8 

7 

32.8 

48 

32.8 

63 

32.8 

5 

32.8 

6.3 

16 

4.8 

114 

•  8 

i  6 

»  3 

i  36 

— 

i  21 

— 

»  10 

— 

i  22 

— 

•  112 

— - 

*  Average  of  10  most  quiet  days. 
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Theodore  R.  Ballard 
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Pergaaon  Press,  Inc. 
Maxwell  House 
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Elmsford,  New  York  10523 
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in  the  Spring  quarter  of  1991. 
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and  attach  any  forms  or  instructions  required  to  confirm  permission.  If  a 
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this  material. 
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